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1.  ABSTRACT 


A  comprehensive  experimental  study  of  the  near  millimeter  wave  dynamics  of  a 
variety  of  solids  is  being  carried  out.  Since  in  many  cases  the  investigation  involves 
driving  the  system  far  from  equilibrium  with  IR  or  UV  laser  radiation,  knowledge  of  the 
time  dependent  optical  response  over  a  broad  frequency  region  is  required.  Effects  such  as 
persistent  IR  spectral  hole  burning  and  population  inversion  in  internal  vibrational  modes  in 
solids  have  been  uncovered  using  this  approach. 

The  discovery  that  persistent  IR  hole  burning  can  be  used  to  manipulate  the  near 
millimeter  wave  response  of  defects  in  crystals,  such  as  N02‘:KI,  and  that  it  also  provides 
a  sensitive  probe  of  the  optical  photodarkening  in  chalcogenide  glasses  have  opened  up 
new  research  avenues  which  require  simultaneous  probes  of  very  different  frequencies  in 
the  IR  and  visible. 

Our  on  going  exploration  of  the  far  IR  response  of  high  temperature 
superconductors  in  sintered,  thin  film,  single  crystal  and  even  single  grain  form  continue  to 
provide  new  information  about  the  dynamical  response  of  these  unusual  materials  both  in 
the  normal  and  superconducting  states. 
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2.  EXECUTIVE  SUMMARY 


Progress  in  the  understanding  of  the  optical  pumping  of  defect  modes,  in  particular, 
and  the  near  millimeter  wave  (NMMW)  properties  of  solids,  in  general,  has  been  pursued 
over  a  wide  front  during  the  past  three  years.  One  measure  of  the  effort  is  the  27  published 
papers,  26  published  talks  and  4  theses  which  have  been  completed. 

During  this  period  the  new  field  of  high  temperature  superconductivity  was 
discovered.  Since  the  energy  gap  for  these  materials  occurs  in  the  middle  of  the  far  infrared 
(FIR)  and  since  the  electromagnetic  skin  depth  is  ~  2000A,  this  technique  provides  an 
important  window  through  which  to  observe  the  dynamics  of  these  novel  mater  als.  The 
other  possible  technique,  electron  tunneling,  generates  data  which  has  proven  to  be 
extremely  difficult  to  interpret  because  the  electrons  probe  the  contaminated  surface  which 
is  not  characteristic  of  the  bulk  since  the  range  of  the  superconducting  wave  function  is 
comparable  to  the  contamination  thickness  ~  10  to  30  A. 

From  our  measurements  on  sintered  material,  single  crystals  and  thin  films  we  have 
been  able  to  show  that  the  intrinsic  properties  of  the  high  Tc  materials  are  still  to  be 
determined.  Still  the  general  behavior  of  the  electromagnetic  response  is  understood 
qualitatively.  These  materials  in  the  superconducting  state  do  have  extremely  high 
reflectivities  in  the  NMMW  region  but  in  the  IR  region  the  reflectivity  is  quite  small  and 
temperature  dependent  similar  to  that  of  nonsuperconducting  ceramic  oxides. 

Our  most  recent  measurements  have  been  on  isolated  single  crystal  grains.  The 
results  show  that  in  the  superconducting  state  a  giant  electric  dipole  resonance  occurs  in  the 
NMMW  region  for  La1.xSrxCu04.y.  The  resonant  frequency  can  be  tuned  throughout  this 
spectral  region  by  simply  changing  x.  These  results  specify  fairly  precisely  the  physical 
parameters  characterizing  the  response  of  the  system  in  both  the  normal  and  the 
superconducting  states. 

A  detailed  discussion  of  our  investigation  into  the  electromagnetic  properties  of 
these  superconductors  is  presented  in  the  Progress  Report  in  Section  (3.1).  Since  the 
results  do  confirm  that  FIR  spectroscopy  is  an  important  tool,  it  is  proposed  that  this  work 
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be  continued  Continued  comparison  of  theory  with  experiment  for  the  sintered  material  is 
extremely  important  because  at  present  such  theoretical  analysis  does  not  give  results 
consistent  with  the  thin  film  and  single  crystal  data.  Since  all  new  high  Tc  superconductors 
first  appear  as  sintered  pellets  this  disagreement  needs  to  be  understood  if  the  FIR 
measuements  are  to  reach  their  full  potential. 

Persistent  IR  spectral  hole  (PIRSH)  burning  has  now  been  observed  for  electronic 
defect  states  in  disordered  semiconductors,  molecules  in  dielectric  crystals  and  molecules  in 
semiconducting  glasses.  We  discovered  that  the  excited  electronic  defect  states  in  ordered 
semiconductors  overlapped  each  other  producing  homogeneously  broadened  transitions 
which  could  not  be  hole  burned.  Disordered  semiconductors,  on  the  other  hand,  had 
inhomogeneously  broadened  transitions  which  did  indeed  show  persistent  effects.  Since 
all  internal  vibrational  modes  of  molecular  defects  in  crystals  are  extremely  well  localized 
the  transitions  are  inhomogeneously  broadened  hence  such  systems  are  good  candidates  for 
high  resolution  PIRSH  burning.  The  effect  has  been  found  for  a  variety  of  symmetry  types 
ranging  from  spherical  to  diatomic  to  asymmetric. 

The  N02':  alkali  halide  system  has  proved  to  be  particularly  rewarding  since  a 
variety  of  new  persistent  effects  in  the  IR  and  the  FIR  have  been  discovered.  Burning  in 
the  IR  produced  the  following  results:  first,  hole  burning  one  part  of  a  tunneling  manifold 
in  the  vibrational  transition  produced  holes  in  other  parts,  second,  hole  burning  one 
vibrational  mode  changes  the  absorption  strengths  of  the  other  internal  vibrational  modes 
and  third,  the  same  procedure  changes  the  FIR  spectrum  of  the  crystal.  This  last  effect  is 
quite  interesting  as  well  as  useful  since  it  demonstrates  that  the  properties  of  the  localized 
phonon  modes  which  involve  both  the  lattice  and  the  molecule  are  directly  coupled  to  the 
internal  modes  of  the  molecule  alone.  Finally  the  internal  vibrational  modes  are  also 
directly  coupled  to  the  local  axes  of  the  molecule  in  the  lattice.  This  coupling  causes  the 
molecule  to  reorient  while  irradiated  with  weak  diode  laser  radiation  at  the  vibrational 
resonance  frequency.  Once  reoriented  the  molecule  must  overcome  a  potential  barrier  to 
return  to  the  equilibrium  configuration.  Since  this  is  an  activated  process,  at  low 
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temperatures  the  recovery  is  suppressed  hence  the  appearance  of  persistent  holebuming  in 
this  system.  This  technique  has  been  used  successfully  to  determine  the  equilibrium 
orientation  of  the  molecule  with  respect  to  the  crystal  axes  of  the  host  lattice. 

The  first  attempt  to  hole  bum  this  defect:lattice  system  directly  in  the  FIR  has  been 
made  possible  by  first  identifying  all  of  the  isotopes  of  the  gap  modes  and  then  fine  tuning 
the  mode  frequency  by  alloying  another  alkali  halide  with  the  host  lattice.  The  mode  can  be 
shifted  to  higher  or  lower  frequencies  by  0.2  cm-1  using  this  technique  of  "poor  man 
tuning".  No  persistent  hole  burning  has  been  observed.  The  negative  result  indicates  that 
either  the  FIR  laser  did  not  have  sufficient  intensity  to  determine  the  low  burning  efficiency 
or  that  the  photon  energy  is  too  small  to  put  the  molecule  in  the  nearly  forbidden  excited 
state  configuration.  Most  of  Section  (3.2)  is  directed  towards  explaining  these 
measurements.  More  work  needs  to  be  carried  out  on  FIR  holebuming  to  determine  which 
limit  is  correct.  In  addition,  our  discovery  that  PIRSH  burning  occurs  in  combination 
bands  in  the  near  IR  also  needs  to  be  examined  more  systematically. 

The  discovery  of  PIRSH  burning  at  molecules  in  chalcogenide  glasses  has  opened 
up  a  completely  different  kind  of  persistence  mechanism.  Now  one  does  not  rely  on  the 
molecule  to  produce  more  than  one  ground  state  configuration  to  support  persistent 
dynamics  but  instead  rearrangement  of  the  disordered  lattice  itself.  The  fact  that  we  have 
observed  persistence  for  molecules  which  form  a  chemical  bond  (HS:  inhomogeneous 
width  ~  200  cm1)  as  well  as  molecules  which  are  simply  trapped  in  the  glass  upon  cooling 
(C02:  ~  5  cm1)  support  this  statememt.  A  description  of  this  work  is  given  in  Section 
(3.2.4).  It  should  be  noted  here  that  the  HS  molecule  produces  a  distinct  absorption  feature 
in  commercial  meter  long  IR  fibers  made  of  AS2S3.  With  sufficient  intensity  it  should  be 
possible  carry  out  PIRSH  burning  in  the  fiber  at  room  temperature  even  though  the 
persistence  time  would  be  extremely  small. 

A  new  research  direction  involves  PIRSH's  in  the  chalcogenide  glasses  as  sensitive 
detectors  of  atom  rearrangemehts.  The  investigation  involves  using  the  interaction  between 
the  photodarkening  (PD)  of  thfese  glasses  by  band  gap  radiation  and  PIRSH  generation  to 
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understand  PD  itself.  Because  PIRSH  burning  centers  are  by  their  very  nature  sensitive  to 
changes  in  their  local  environment  any  process  which  produces  long  range  structural 
changes  in  the  glassy  host  should  produce  measurable  changes  in  the  persistent  holes. 
Recently,  we  have  searched  for  such  an  effect ,  and  found  that  extremely  large  changes  in 
holes  can  be  produced  with  surprisingly  low  doses  of  band  gap  radiation.  We  intend  to  set 
up  an  interference  pattern  with  a  band  gap  laser  and  study  persistent  hole  production  as  a 
function  of  fringe  spacing  to  investigate  the  range  of  influence  of  the  photo-induced 
changes.  The  results  would  help  decide  whether  PD  is  associated  with  migrating  localized 
defects  or  the  long  range  rearrangement  of  many  atoms. 

In  1988  it  was  suggested  by  the  PI  and  his  colaborator  Professor  S.  Takeno  in 
Japan  that  plane  wave  phonons  were  insufficient  to  describe  the  thermodynamic  properties 
of  crystals  at  finite  temperature.  It  was  proposed  that  in  addition  to  the  plane  wave  modes 
which  minimized  the  lattice  anharmonicity  the  crystal  also  contained  localized  vibrational 
modes  which  maximized  the  crystal  anharmonicity.  These  localized  modes  move  as 
particles  which  have  an  internal  vibrational  degree  of  freedom.  We  call  them  "anons". 

The  new  feature  with  respect  to  the  thermodynamic  properties  of  glasses  is  that  they 
diffuse  through  the  lattice  and  do  not  transmit  sound  like  phonons  do.  We  have  shown  that 
these  modes  can  account  for  the  unusual  specific  heat  of  glasses  since  the  translational 
diffusive  motion  produces  a  linear  temperature  dependent  contribution  to  the  specfic  heat 
and  the  internal  vibrational  degree  of  freedom  provides  the  extra  cubic  temperature 
dependence,  over  an  above  Debye,  both  of  which  are  observed  in  experiments.  A  more 
detailed  description  of  these  ideas  is  presented  in  Section  (3.3). 

Although  we  have  shown  that  the  presence  of  anons  can  describe  the  low 
temperature  specific  heat  of  glasses,  we  propose  that  these  particles  could  exist  in  doped 
crystals  as  well.  The  localized  strain  field  near  an  extrinsic  defect  (such  as  an  ion  or  a 
molecule)  could  act  as  an  anon  trap.  We  are  investigating  the  possibility  that  the  Ag+ 
impurity  in  KI:Ag+  may  have  such  an  anon  bound  to  it  at  low  temperatures  but  that  the 
excitation  may  be  released  as  the  temperature  is  increased  to  20  K.  This  process  could 
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explain  the  n^id  transformation  of  the  Ag+  ion  from  an  on-center  to  an  off-center  impurity 
as  the  temperature  is  varuied  over  this  limited  range. 

Another  place  where  a  trapped  anon  may  influence  experiment  is  in  the  incoherent 
saturation  of  the  internal  vibrational  modes  of  molecules  in  crystals.  Exciting  the  internal 
mode  of  the  molecule  could  provide  enough  local  energy  to  release  the  anon  from  the  trap 
site.  This  extra  degree  of  freedom  would  provide  another  relaxation  channel  besides 
photon  emission  or  multiphonon  emission  by  which  a  vibrational  quantum  could  be 
removed  from  the  molecule. 

The  arrival  of  a  tunable  UV-IR  Spectra  Physics  laser  gives  us  the  unique 
opportunity  to  pump  the  electronic  and  vibrational  transitions  of  different  point  defects  and 
molecules  while  at  the  same  time  monitoring  their  translational  vibrational  spectrum  in  the 
NMMW  region.  One  expected  effect  at  low  temperatures  is  persistent  UV  hole  burning. 
The  induced  change  is  expected  to  alter  the  impurity  induced  absorption  spectrum  in  the  far 
infrared.  In  addition,  this  pump  will  permit  the  localized  and  resonant  vibrational  modes 
associated  with  the  defect  to  be  driven  in  a  controlled  manner  far  from  equilibrium. 
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3.  PROGRESS  REPORT 

3.1  The  Electrodynamics  of  High  Temperature  Superconductors 
3.1.1  Background  information 

The  far-IR  properties  of  the  high  Tc  compounds!  1]  in  sintered  polycrystalline,  thin 
film  and  single  crystal  form  continue  to  surprise  and  fascinate  solid  state  spectroscopists  in 
ever  increasing  numbers.  Although  the  experimental  situation  is  in  constant  turmoil  nearly 
everyone  agrees  that  in  the  visible  the  material  is  black  as  coal  with  a  low  reflectivity  like 
that  of  a  wide  band  gap  semiconductor  while  in  the  far  infrared  the  reflectivity  in  the 
superconducting  state  corresponds  to  that  of  a  metal  not  a  lossless  metal,  characteristic  of  a 
superconductor  at  frequencies  below  the  energy  gap,  but  still  a  good  metal.  Divining  more 
than  this  from  the  data  is  what  continues  to  produce  controversy  but,  on  the  other  hand,  it 
is  also  what  makes  these  materials  of  spectroscopic  interest  To  appreciate  the  fundamental 
IR  problems  it  is  useful  first  to  review  the  electromagnetic  response  for  low  temperature 
superconductors. 

In  general,  a  nonlocal  relation  between  the  current  density  and  the  electric  field 
together  with  Maxwell's  equations  must  be  used  to  describe  the  fields  inside  a  pure  metal  at 
low  temperatures.  The  surface  impedance  concept  permits  a  reformulation  of  the  problem 
in  terms  of  a  boundary  condition  on  the  fields  at  the  metal  surface.  The  defining  equation  is 

Et  =  z  Htx  n  (1) 

where  the  subscript  t  signifies  the  external  tangential  fields  evaluated  at  the  boundary,  the 
normalized  surface  impedence  z  =  7/Zq  with  Zq  the  impedence  of  free  space  and  n  is  the 
normal  to  the  metal  surface.  The  complex  coefficient  z  completely  characterizes  the 
medium  so  the  normal  incident  reflectivity  R  and  the  absorptivity  A  are  related  by 

A  =  1  -R  =  - (2) 

(1  +  r)2  +  x2 

For  most  metals  r,  x  «  1  so  that  the  absorptivity  reduces  to 

A  =  4r.  (3) 

The  absorptivity  of  a  dirty  metal  in  the  far  infrared  corresponds  to  the  classical  skin 
effect  limit.  In  this  regime  the  electromagnetic  skin  depth  8  is  large  compared  to  an  electron 
mean  free  path  vft  so  that  a  local  conductivity  o(co)  can  be  defined.  If,  further,  it  is 
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assumed  that  that  the  carrier  scattering  time  is  independent  of  energy  and  that  the  Drude 
model  can  be  used  to  describe  the  a.c.  conductivity  then 

0(G))  =  0O(1  -icot)-1  ,  0Q  =  (C0p2  X)  /  47C  (4) 

where  0O  is  the  dc  conductivity  and  cop  is  the  plasma  frequency.  Note  that  two  parameters 
x  and  G)p  characterize  the  electromagnetic  properties  of  the  metal.  The  corresponding 
components  of  the  surface  impedance  are 

r  =  (l/2oopx)(2cDX)1/2[(l  +  co2!2)1^  -  cut]1/2  (5) 

and 

x  =  -  (l/2copx)(2cot)1/2[(l  +  co^c2)1/2  +  cox] 1/2  .  (6) 

In  the  far  infrared,  cox  «  1  and  the  absorptivity  at  normal  incidence  has  the  simple 
Hagen-Rubens  form 

Ah(co)  =  (2co  /  7t  0O)1/2  (7) 

in  which  An  increases  as  to1/2.  At  a  frequency  of  100  cm1  and  a  resistivity  of  70  (jXicm, 
An  =  0.03.  For  larger  infrared  frequencies,  in  the  limit  where  1  «  cox  «  copx,  the 
absorptivity  has  the  constant  value 

An(co)  =  2  /  (cop  x)  .  (8) 

At  still  larger  frequencies,  in  the  visible,  transitions  to  other  bands  become  important  and 
the  simple  intraband  absorption  model  given  here  is  no  longer  applicable. 

For  materials  with  small  conductivities  the  high  frequency  interband  transitions  can 
still  influence  the  the  size  of  the  absorptivity  in  the  IR  and  FIR  region  through  their  effect 
on  the  displacement  current  in  the  medium.  This  effect  can  be  seen  most  easily  when  the 
dielectric  function  of  the  medium  e(co)  is  written  in  terms  of  the  contribution  from  the 
displacement  and  conduction  currents  so 

e((o)  =  £„  +  4rci  o(co)/co  (9) 

where  ex  is  the  constant  low  frequency  contribution  from  the  interband  transitions.  When 
this  contribution  is  significant  it  is  often  useful  to  plot  0j(co)  and  e^G))  rather  than  0j(g)) 
and  02(co)  so  that  the  two  Drude  parameters  x  and  G>p  will  be  evident  If  the  material  has  a 
magnetic  response  as  well  as  an  electric  one  then  to  compare  theory  and  experiment  in  the 
local  limit  the  connection  between  the  surface  impedance  and  the  different  susceptibilities 
is 
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z  =  [|j.(co)/e(co)] 1/2  .  (10) 

The  dotted  curve  in  Fig.  (1)  shows  the  measured  surface  resistance  at  room 
temperature  for  an  annealed  lead  foil  in  the  far  infrared  spectral  region[2].  A  multiple 
parallel  plate  waveguide  technique  was  employed  to  measure  this  small  absorptivity.  In 
this  frequency  region  the  displacement  current  is  negligible  compared  to  the  conduction 
term  and  also  p,  =  1.  The  solid  curve  in  the  figure  is  calculated  from  the  Drude  theory  with 
the  measured  dc  conductivity  of  the  foil  inserted  into  Eq.  (7).  With  no  free  parameters 
good  agreement  is  obtained  between  theory  and  experiment  throughout  this  spectral  region. 

Because  of  Landau  damping  the  FIR  absorptivity  of  a  pure  metal  at  low 
temperatures  no  longer  agrees  with  the  simple  Drude  predictions  based  on  a  temperature 
dependent  relaxation  time.  This  anomalous  skin  effect  regime  has  been  well  studied  both 


Figure  1. 


in  the  normal  state  and  in  the  superconducting  state.  The  frequency  dependent  response  of 
a  superconductor  in  the  extreme  anomalous  limit  as  calculated  by  Mattis  and  Bardeen[3]  is 
shown  in  Fig.  (2).  The  electromagnetic  gap  fuog  =  2A  =  3.5  kBT  for  this  weak  coupling 
limit.  The  normal  state  conductivity  an  is  assumed  to  be  entirely  real  in  this  frequency 
interval  (cm  «  1)  and  the  results  are  shown  as  a  ratio  of  the  real  at  and  imaginary  o2  in  the 
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superconducting  state  to  cn.  These  conductivity  expressions  also  apply  in  the  opposite  or 
extreme  dirty  limit.  The  ratio  of  the  extreme  anomalous  surface  resistance  in  the 
superconducting  state  to  the  normal  state  for  a  bulk  sample  is  also  shown  in  Fig.  (2).  For 
the  opposite  limit  Tinkham[4]  has  shown  that  the  frequency  dependence  of  the  surface 
resistance  remains  qualitatively  the  same  as  in  Fig.  (2). 


For  a  strongly  coupled  superconductor  such  as  Pb  the  spectrum  for  the 
superconducting  state  is  made  complex  by  the  appearence  of  the  Holstein  absorption 
process[5-8].  The  coupling  between  electrons  and  phonons  in  metals  makes  allowed  an 
infrared  absorption  process  at  low  temperatures  associated  with  phonon  generation  by  the 
excited  electrons.  The  product  of  the  square  of  the  electron-phonon  matrix  element  and  the 
phonon  density  of  states  in  the  metal  can  be  obtained  from  a  measurement  of  the  difference 
in  the  FIR  absorptivity  between  the  superconducting  and  normal  state. 

Measurements^]  of  the  difference  in  the  FIR  surface  resistance  between  super  and 
normal  for  the  extreme  anomalous  (Pb)  and  the  extreme  dirty  (Pb-O.lln)  limits  are  shown 
in  Fig.  (3).  Above  the  gap  frequency  the  absorptivity  in  the  superconducting  state  exceeds 
that  in  the  normal  state  up  to  frequencies  about  10  x  the  gap  frequency.  Note  that  the  alloy 
sample  has  a  much  larger  surface  resistance  in  the  normal  state  below  the  gap  than  does  the 
pure  Pb  sample  and  that  it  also  has  a  corresponding  larger  surface  resistance  in  the 
superconducting  state  above  the  gap.  These  results  illustrate  the  "BrSndli  sum  rule"  that  the 
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absorptivity  in  the  normal  state  below  the  gap  must  equal  the  excess  absorptivity  in  the 
superconducting  state  above  the  gap[9],  i.e.. 


Another  interesting  feature  about  the  data  shown  in  Fig.  (3)  is  that  the  frequency 
dependence  of  the  excess  absorptivity  in  the  superconducting  state  is  completely  different 
for  the  two  cases.  For  pure  Pb  the  profile  is  associated  with  the  availablility  of  phonons 
hence  for  T  =  0  with  the  spectrum  of  the  phonon  density  of  states  whereas  for  the  alloy 
case  the  main  ingredient  is  quasi-elastic  scattering  near  zero  frequency. 


3.1.2  High  Tc  sintered  Samples 

The  room  temperature  FIR  reflectivity  of  sintered  polycrystalline  YBa2Cu307.y  >s 
represented  by  the  solid  line  in  Fig.  (4a)[10].  The  classical  skin  effect  prediction  (dashed 
line)for  a  dc  resistivity  of  800  p£2cm  (the  measured  resistivity  is  890  (iflcrn  ±  10%)  is  also 
shown.  At  the  lowest  frequency  limit  of  our  data  the  two  curves  converge,  but  for  higher 
frequencies,  the  reflectivity  deviates  strongly  from  the  Drude  prediction.  We  also  observe 
very  little  temperature  dependence  to  the  reflectivity  ( 100  K,  dotted  line)  in  agreement  with 
other  measurements  at  larger  frequencies.  While  the  dc  resistivity  changes  by  a  factor  of 
three  from  just  above  Tc  to  room  temperature,  the  reflectivity  changes  by  only  a  few 


11 


percent[10]. 

Figure  (4b)  shows  the  room  temperature  reflectivity  as  a  function  of  frequency  on  a 
logarithmic  frequency  axis,  covering  four  decades  in  frequency.  These  results  are  similar 
to  those  reported  by  others[l  1-13].  The  appearance  of  phonon  modes  such  as  those  shown 
in  the  100  to  1000  cm*1  region  is  not  common  in  metallic  reflectivity,  but  they  have  been 

seen  in  materials  with  small  optical  conductivity!  14]. 
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How  can  one  analyze  such  sintered  material  data  more  quantitatively?  With  the 
definitive  observation!  15]  in  the  FIR  reflectivity  of  2-dimensional  optical  conductivity  in 
the  ab  plane  of  single  crystal  La2Ni04,  a  compound  structurally  isomorphic  to 
Lai .ssSro.isCuO^y  (the  40  K  superconductor),  an  effective  medium  modeling  of  the  FIR 
data  for  both  the  superconducting  and  normal  state  of  the  sintered  material  becomes 
practical.  For  the  grain  size  small  compared  to  the  wavelength  and  for  2-D  conduction  with 
the  third  direction  in  the  grain  insulating,  the  effective  medium  approximation  (EM A) 
equation  has  a  relatively  simple  form,  with  a  fixed  fill  fraction  of  f  =  2/3  for  the  metallic 
component.  To  make  the  problem  tractable  we  have  assumed  that  the  crystallites  are 
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spheroids  with  the  same  symmetry  as  the  conductivity  tensor,  so  that  a  single  parameter, 
the  depolarization  factor,  L,  along  the  c  direction  of  the  crystallite,  specifies  the  EMA.  The 
effective  conductivity  for  the  medium  is  obtained  by  solving  the  following  equation[16]: 

1 _ - +  1 - °C~CT! - =  0  (12) 

3  (l+L)oe+  (l-Do^b  3  2[(l-L)oe  +  Loc] 

where  oe  is  the  desired  conductivity  of  the  isotropic  effective  medium  in  terms  of  the 
known  conductivities  of  the  individual  anisotropic  grains. 

For  the  case  of  spherical  crystallites,  L  =  1/3  and  Eq.  (12)  reduces  to  the  standard 
EMA  result  for  an  anisotropic  polycrystalline  material.  For  needle-like  crystallites  L  =  0 
along  the  c  axis,  while  for  plate- like  crystallites  L  =  1.  We  find  that  this  EMA  with  Lorentz 
oscillator  contributions  for  the  phonons  and  a  two  dimensional  Drude  free  carrier 
component  provides  a  good  description  of  the  normal  metal  results  for  the  sintered  systems 
throughout  the  FIR  region[16,17]. 

Figure  5(a)  presents  our  FIR  data  on  the  reflectivity  of  Laj  85Sr0 15Cu04.y  in  the 
normal  state  (T  =  40  K,  dashed  curve)  and  in  the  superconducting  state  (T=10  K,  solid 
curve)[10].  The  dashed  line  fit  shown  in  Fig.  5(b)  demonstrates  that  the  Drude  term  does 
not  provide  enough  reflectivity  at  low  frequencies.  The  calculated  reflectivity  for  the 
superconducting  state  represented  by  the  solid  line  in  Fig.  5(b)  makes  use  of  the  Mattis- 
Bardeen  equations  to  model  the  2D  superconductivity  in  the  a,b  plane.  This  too  reproduces 
some  of  the  qualitative  features  observed  in  Fig.  5(a),  for  example,  although  2D  is  chosen 
as  68  cm1,  the  onset  of  absorption  occurs  at  about  35  cm1.  Also  the  reflectivities  cross  at 
73.7  cm1  and  the  normal  state  becomes  more  reflecting  than  the  superconducting  state 
above  this  frequency. 

To  demonstrate  the  key  role  of  grain  geometry,  Fig.  5(c)  shows  the  reflectivity  in 
the  superconducting  and  normal  states  as  calculated  with  the  EMA  using  identical 
parameters  as  Fig.  5(b),  except  that  Lc  =  0.7,  corresponding  to  oblate  spheroids. 
Immediately  obvious  is  the  dramatic  loss  of  the  enhanced  absorptivity  in  the 
superconducting  state  above  the  gap  and  the  shift  in  position  of  the  equal  reflectivity 
frequency  which  often  has  been  associated  with  the  frequency  of  the  energy  gap,  2A. 
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Figure  5. 


1.0 


R  0.8  h  V 


0.6 


(c) 


0  100  200  300 

Frequency  (cm-1) 


The  qualitative  agreement  between  the  EM  A  calculations  and  the  measured 
reflectivity  on  the  sintered  material  in  the  long  wavelength  limit  demonstrates  that  the 
individual  grains  are  highly  anisotropic  in  the  far  infrared  spectral  region.  The  good 
specular  properties  of  the  sintered  samples  in  the  visible  indicate  that  the  grains  must  be 
nearly  isotropic  at  these  wavelengths.  Between  these  two  limits  the  grain  size  is 
comparable  to  the  wavelength  of  the  radiation  and  the  grains  are  anisotropic  so  that  the 
surface  appears  rough  and  scatters  radiation.  So  far  there  is  no  theory  to  describe  the 
results  in  this  region. 

3.1.3  Optical  Properties  of  composites 

Recently  there  has  been  much  controversy  surrounding  the  origin  and  significance 
of  a  peak  centered  around  0.5  eV  in  the  optical  conductivity  of  sintered  high  temperature 
superconductors.  Several  reports[18,19]  have  been  published  on  the  IR  spectrum  of 
sintered  YBa2Cu307_x  and  La2.xSrxCu04  which  claim  that  the  observed  correlation 
between  the  oscillator  strength  of  the  "0.5  eV"  peak  and  the  degree  of  "superconductivity" 
in  the  samples  might  indicate  that  pairing  in  these  materials  is  mediated  by  a  high-energy 
r'^tronic  excitation  such  as  an  exciton[20].  Later,  IR  measurements  of  single  crystals 
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were  made  by  many  workers,  using  both  reflectivity[21,22]  and  direct  absorption 
methods[23].  Most  of  these  measurements  on  single  crystals  have  failed  to  show  a  strong 
"0.5  eV"  peak  although  a  broad  absorption  is  observed  to  extend  throughout  the  IR  spectral 
region.  Figure  (6)  shows  an  IR  peak  in  the  conductivity  for  three  different  sintered  copper 
oxide  compounds.  The  top  trace  is  for  the  40  K  superconductor,  the  middle  trace  is  for  the 
90  K  superconductor  and  the  bottom  trace  is  for  a  non-superconductor[24].  The 
compound  La4BaCu5013  can  be  viewed  as  a  cubic  perovskite  with  an  array  of  channels  of 
oxygen  vacancies  running  parallel  to  the  c  axis,  so  the  Cu-0  network  is  three  dimensional 
but  anisotropic [25].  The  3-D  conductivity  screens  the  phonon  dipole  moments  and  keeps 
the  lattice  component  from  showing  up  even  in  the  sintered  pellets. 


Figure  6. 
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Why  the  "0.5  eV  peak"  appears  in  the  sintered  copper-oxide  perovskite-like 
materials  is  an  interesting  question  even  if  it  is  not  an  intrinsic  effect.  Possible  answers 
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could  involve:  (1)  effects  of  the  composite  nature  of  the  system  such  as  sphere  resonances; 
or  (2)  the  strong  anisotropic  nature  of  the  material  and  the  random  orientation  of  the  grains. 
In  a  recent  Comment,  Orenstein  et  al., [26]  proposed  a  simple  model  which  explains  the 
"0.5  eV  absorption  peak"  in  terms  of  the  optical  anisotropy  of  the  material  in  the 
geometrical  optics  limit  While  this  model  seems  to  reproduce  qualitatively  the  presence  of 
the  peak  in  the  mid-IR,  we  believe  that  the  complete  picture  must  be  more  complicated.  In 
this  region  of  the  spectrum,  the  size  of  the  anisotropic  grains  is  such  that  neither  the 
geometrical  nor  effective  medium  limits  can  be  correct,  and  scattering  effects  must  be 
important.  It  may  be  that  the  origin  of  the  peak  can  be  determined  more  clearly  by  careful 
studies  in  which  the  grain  size  and  porosity  of  the  materials  can  be  controlled. 

3.1.4  Films  and  Single  Crystal  Oxide  Superconductors 

It  is  useful  to  compare  the  superconducting  response  of  sintered,  thin  film  and 
single  xtal  FIR  materials.  Shown  in  Figure  7(a)  is  the  ratio  of  the  reflectivity  of  sintered 
YBa2Cu307.y  in  the  superconducting  state  (10  K)  to  that  in  the  normal  state  (100  K).  The 
polycrystalline  nature  is  accounted  for  again  with  the  EMA[27].  The  dashed  curve  in  Fig. 
7(a)  is  obtained  from  such  a  fit  to  the  normal  state  data,  followed  by  the  application  of  the 
MB  equations  to  model  the  free  electron  conductivity.  Here  a  gap  of  170  cm1  was 
obtained  for  plate-like  crystallites.  An  equally  satisfactory  fit  is  obtained  with  a  gap  of  150 
cm1  for  spherical  crystallites.  The  relative  superconducting  gap  energies  obtained  from 
sintered  sampes  of  both  YBa2Cu307.y  and  ^Sro.^CuO^y  are  thus  seen  to  be  similar, 
with  2A/kTc  =  2.6. 

Contrasting  with  the  behavior  of  the  bulk  samples.  Figure  7(b)  shows  data  for  1 
pm  thick  films  of  YBa2Cu3(>7.y  on  Zri)2.  These  films  are  of  high  quality,  having  Tc  =  85 
K  and  Jc  >  106  A/cm2  at  T=  4.2  K.  The  films  are  somewhat  oriented,  with  the  c-axis 
preferentially  perpendicular  to  the  film  surface.  Using  the  same  dielectric  function  in  the  a- 
b  plane  obtained  from  fitting  the  data  on  bulk  material,  but  with  the  measured  dc  resistivity 
of  200  pficm  in  the  a-b  plane,  produces  the  dashed  curve  in  Fig.  7(b).  Here  the  choice  of 
the  gap  energy  indicates  2A/kTc  =  6.4  ,  larger  than  the  value  obtained  for  bulk  sintered 
material.  The  conductivity  obtained  from  a  Kramers-Kronig  transformation  of  the 
reflectivity  covering  over  three  orders  of  magnitude  in  frequency  does  not  show  the 
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electronic  mode  in  the  IR  found  for  bulk  polycry  stalling  samples. 


Figure  7. 


Just  as  the  relative  gap  parameter  for  the  films  exceeds  that  of  the  sintered  material, 
so  that  of  the  single  crystals  exceeds  that  of  the  films.  Single  crystals  have  been  pieced 
together  to  form  a  mosaic  with  the  c-axis  normal  to  the  surface.  Our  reflectivity 
measurements  shown  in  Fig.  7(c)  and  our  determination  of  the  gap  parameter,  2A/kTc  =  8, 
are  in  agreement  with  those  reported  on  single  crystals  by  others  [21].  Again  the  dashed 
curve  is  a  fit  to  the  data  using  the  dielectric  function  determined  on  the  bulk,  with  a  dc 
resistivity  of  200  pQcm.  The  gap  value  is  determined  from  this  fit.  As  for  the  films,  the 
reflectivity  signature  of  the  superconducting  state  is  not  "complete"  in  that  there  is 
absorption  of  energy  below  2A.  Thus  the  dashed  curve  in  Fig.  7(c)  has  been  reduced  by 
50%  to  agree  with  the  overall  magnitude  of  the  data.  Our  measurements  on  different  single 
crystals  show  that  this  absorption  below  the  gap  is  not  the  same  from  sample  to  sample 
hence  it  must  be  an  extrinsic  effect. 

Recently  Timusk  et  a/., [28]  proposed  an  alternate  explanation  for  the  Rj/R,,  results 
which  does  not  make  specific  use  of  a  gap  or  equivalently  assumes  that  X'1  «  2A  so  that 
the  entire  free  carrier  Oj((d)  profile  in  the  normal  state  is  compressed  into  the  delta  function 
at  ©  =  0  in  the  superconducting  state.  With  x1  small  and  Oj(dc)  known  for  the  normal 
state,  the  conductivity  sum  rule  gives  a  small  plasma  frequency.  From  Hall 
measurements[29]  an  estimate  of  the  density  of  free  carriers  can  be  obtained  independently. 


17 


This  value  together  with  the  small  plasma  frequency  implies  a  large  carrier  mass;  however, 
specific  heat  gives  a  carrier  mass  which  is  not  particularly  heavy [30].  At  present  all  of 
these  results  do  not  give  a  consistent  picture. 

Currently  the  far  infrared  (FIR)  spectroscopic  data  on  high  Tc  single  crystals  in  the 
superconducting  state  can  be  interpreted  as  evidence  for  a  large  energy  gap[31],  a  BCS-like 
gap[32]  or  no  gap  at  all[33].  The  optical  properties  of  crystals  of  variable  quality  make  it 
difficult  to  conclude  whether  or  not  the  peak  observed  in  the  ratio  of  the  FIR  reflectivity  of 
the  ab  plane  of  YBa2Cu307.y  in  the  superconducting  state  to  that  of  the  normal  state  should 
be  interpreted  as  evidence  for  the  energy  gap[31,32]  or  as  evidence  for  a  zero  crossing  of 
the  real  part  of  the  dielectric  function[33].  The  gap  of  La2_xSrxCu04.y  is  even  less  well 
determined  because  it  has  been  difficult  to  grow  single  crystals  at  the  Sr  doping  which 
produces  the  highest  Tc's.  The  optical  properties  of  superconducting  sintered  La2_ 
xSrxCu04_y  have  been  analyzed  in  some  detail[34]  but  the  excess  absorption  in  the 
superconducting  state  associated  with  transport  between  grains  obscures  the  determination 
of  a  true  gap  signature.  Compared  to  the  extensive  work  on  these  polycrystalline[34,35]  or 
single  crystal  samples[31-33],  the  FIR  electrodynamics  of  high  Tc  superconductors  in  the 
isolated  particle  form  have  been  ignored  even  though  the  relation  between  the  bulk  dielectric 
function  and  the  small  particle  response  is  well  known  and  has  been  investigated  for  many 
materials[36], 

3.1.5  Isolated  Superconducting  Particles 

Our  most  recent  investigation  has  focused  on  the  FIR  response  of  small  isolated 
La2.xSrxCu04.y  particles  in  the  normal  and  superconducting  states[37].  A  far  infrared 
superconducting  sphere  resonance  polarized  along  the  c-axis  has  been  observed  to  decrease 
in  frequency  with  increasing  temperature,  increasing  magnetic  field  or  decreasing  Sr 
doping.  No  corresponding  resonance  is  found  in  the  normal  state.  These  observations 
indicate  that  an  energy  gap  does  exist  in  high  Tc  materials  since  the  energy  gap  frequency 
must  be  larger  than  the  c-axis  sphere  resonance  value.  These  results  will  be  described  by 
Dr.  Noh  in  an  invited  paper  at  the  March"Solid  State"  meeting  of  the  American  Physical 
Society. 

Before  turning  to  the  experimental  data  we  summarize  the  different  kinds  of  spectral 
signatures  which  can  be  expected  to  occur  for  small  conducting  particles.  Let  us  consider 
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two  specific  examples  to  identify  qualitatively  the  spectral  signature  for  superconductors  in 
small  particle  form.  Assume  that  a  BCS  energy  gap  separates  the  superconducting  ground 
state  from  the  single  particle  excited  states,  that  copx  »1  and  that  either  (2A)x  »  1  or 
(2A)x  «  1  for  a  bulk  low  temperature  superconductor.  For  the  first  case  the  entire  normal 
state  conductivity  spectrum  collapses  into  the  delta  function  at  zero  temperature  with 
strength  Aj  =  /  an(co)d(0  =  top2/8.  For  the  second  case,  <Tj(£d)  is  zero  up  to  a  gap 
frequency  2A,  above  which  it  rises  monotonically  to  the  frequency  independent 
conductivity  of  the  normal  state,  an,  since  co  is  still «  x1.  The  strength  An  of  the  8(co) 
term  comes  from  the  missing  area  below  the  gap  frequency [3 8],  i.e.,  An  »  2(2A)on.  We 
show  the  expected  results  for  the  corresponding  sphere  resonances  for  these  two  cases  at 
the  top  of  Table  1  and  below  those  the  results  for  other  possible  sets  of  parameters. 
Depending  on  the  relative  size  of  o)p,  x*1,  and  2A,  various  sphere  resonance  signatures  can 
be  expected  in  small  metallic  particles  in  the  normal  and  superconducting  states. 


Table  1. 


Case 

Conditions 

Observable? 
(normal)  (s-c.) 

Characteristics 

I 

£opt »  1,  (2A)x  >  1 

yes 

yes 

The  line  width  will  narrow  drastically 
upon  cooling  below  Tc. 

II 

£l)pT  »  1,  (2A)t  «  1 

yes 

yes 

The  line  width  will  not  change 
upon  cooling  below  Tc. 

III 

0)pT  £  2.  (2A)t  >  1 

no 

yes 

The  peak  position  will  not  have  a 
strong  temperature  dependence. 

IV 

CDpT  £  2,  (2A)I «  1 
and  eaK  <  2A 

no 

yes 

The  peak  position  will  have  a  strong 
temperature  dependence. 

V 

copx  £  2,  (2A)t  «  1 
and  cow  »  2A 

no 

no 

The  peak  will  be  overdamped  since 
11X1(6(0),,)  +  2eh)  is  large  for  both 
states. 

The  FIR  temperature  dependence  of  the  absorption  coefficient,  a,  for  1% 


Lai  85Sro.i5Cu04.y  *n  Teflon  is  shown  in  Figure  (8a).  The  absorption  coefficient 
difference  Aa(T)  =  [a(T)  -  a(T=40K)]  is  plotted  versus  temperature.  The  low  temperature 
absorption  peak,  located  at  54  cm1,  shifts  to  smaller  frequencies  and  weakens  as  the 
temperature  approaches  Tc.  No  resonance  is  seen  in  the  normal  state.  At  the  resonant 
frequency  the  absorption  in  the  superconducting  state  is  larger  than  the  absorption  in  the 
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normal  state  but  far  from  the  resonance  the  reverse  is  true.  A  six  T  magnetic  field  was 
observed  to  decrease  the  line  strength  by  about  40%  and  produce  effects  similar  to  those 
obtained  with  increasing  temperature.  These  observations,  which  are  consistent  with  Case 
IV  in  Table  1,  indicate  that  the  resonance  is  associated  with  the  appearance  of 
superconductivity. 


Figure  8. 


The  Sr  concentration  dependence  and  hence  the  conductivity  dependence  of  the 
absorption  feature  is  shown  in  Figure  (9).  Since  the  absorption  for  the  same  material 
between  0.3%  and  3%  fill  fraction  scales  with  f,  the  difference  between  the  absorption 
coefficients  in  the  superconducting  and  normal  states  are  normalized  by  the  fill  fraction,  f, 
of  the  superconductor.  The  peak  position  and  the  linewidth  of  the  absorption  line  increase 
monotonically  with  x.  However,  the  normalized  absorption  value  at  the  peak  position  and 
the  strength  of  the  resonant  feature  increase  at  low  concentration  of  Sr  but  reach  a 
maximum  around  x=0.18.  Our  samples  with  x  =0.225  and  0.25  show  no  absorption  lines 
up  to  180  cm*1,  although  they  are  superconducting.  Also,  no  line  is  observed  in  non¬ 
superconducting  x=0.02  and  x=0.04  samples. 

The  general  behavior  of  the  peak  position  and  strength  as  a  function  of  Sr 
concentration  can  be  explained  as  a  change  in  the  value  of  an.  If  on  along  the  c-axis 
increases  with  Sr  doping  like  the  d.c.  conductivity  of  sintered  materials,  the  position  and 
the  strength  of  the  resonance  should  increase.  The  disappearance  of  the  resonant  feature 
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above  x=0.225  probably  results  from  the  resonant  frequency  becoming  much  higher  than 
2D,  so  that  the  limit  of  Case  V  in  Table  1  is  in  effect.  The  lack  of  a  resonant  feature  below 


Figure  9. 


Frequency  (cm-’) 


x=0.04  is  due  simply  to  the  lack  of  superconductivity  in  those  samples. 

In  summary,  we  have  observed  a  unique  resonance  phenomenon  in  well  dispersed 
superconducting  La2_xSrxCu04_y  particles  in  teflon.  The  observed  strength  and  position  of 
the  resonance  identifies  for  the  first  time  the  value  of  the  c-axis  conductivity  for  large  Sr 
concentrations.  The  observed  behavior  of  the  resonant  feature  fixes  the  inequalities  that  the 
normal  free  carrier  and  superconducting  parameters  must  satisfy,  namely,  the  normal  state 
carriers  along  the  c-axis  are  heavily  damped  and  the  c-  axis  energy  gap  frequency  is  larger 
than  the  sphere  resonance  frequency  but  smaller  than  the  normal  state  carrier  scattering  rate. 

All  of  the  data  presented  here  show  that  the  FIR  properties  of  the  high  Tc 
superconductors  do  not  correspond  to  the  low  temperature  ones  described  in  the 
background  section.  The  sintered  Cu-0  materials  can  be  described  fairly  acurately  with  the 
EMA  but  the  sphere  resonance  results  indicate  that  the  derived  gap  value  is  too  small.  At 
higher  frequencies,  the  IR  peak  in  the  effective  optical  conductivity  of  the  sintered  materials 
is  an  interesting  effect  which  occurs  in  a  region  of  the  spectrum  where  composite  media 
analysis  is  not  valid.  The  IR  peak  in  the  conductivity  does  not  appear  to  be  intimately 
connected  to  superconductivity  since  it  appears  in  both  superconducting  and  normal  Cu-0 
materials.  A  detailed  analysis  of  the  single  crystal  data  is  still  complicated  by  the  variability 
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of  the  materials  and  by  the  as  yet  to  be  determined  influence  of  the  heavily  twinned  a,b 
plane  on  the  optical  constants.  In  principle,  the  influence  of  this  single  crystal  twinning  (on 
the  scale  of  1000  A)  on  the  IR  and  FIR  properties  could  be  described  with  still  another 
composite  medium  representation. 
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3.2  Optically  Pumped  Defect  Modes  in  Dielectric  Solids 
3.2.1  Hole  burning  history 

The  formation  of  persistent  spectral  holes  in  electronic  bands  in  solids  was  first 
observed  in  1974  and  at  present  persistent  features  caused  by  both  photochemical  and 
nonphotochemical  processes  continue  to  be  found  in  an  ever  increasing  number  of 
solids[39].  In  the  last  few  years  persistent  IR  hole  burning  has  been  observed  in  the 
vibrational  spectrum  of  molecular  defects  in  crystals,  i.e.,  no  electronic  excitation  is 
involved  and  only  vibrational  degrees  of  freedom  are  photoexcited  with  IR  radiation[40- 
45]. 

Persistent  IR  spectral  holes  provide  a  new  high  resolution  technique  with  which  to 
explore  the  properties  of  defects  in  solids.  The  inhomogeneously  broadened  absorption 
lines  associated  with  the  vibrational  modes  of  matrix  isolated  molecules  are  in  general  not 
coincident  with  the  frequencies  of  ir-line  tunable  gas  laser :  consequently,  the  high  power 
techniques  of  incoherent  saturation  and  transient  hole  burning  cannot  be  applied  to  identify 
the  dominant  relaxation  process  in  most  molecular-defect-lattice  combinations.  With  low 
power  tunable  diodes  two  techniques  are  potentially  available:  vibrational  fluorescence, 
which  gives  the  excited  state  lifetime;  and  persistent  hole  burning  which  not  only  provides  a 
high  resolution  probe  of  the  defect  but  also  identifies  the  excited  state  dephasing  time  of  the 
vibrational  mode. 

What  are  the  most  general  requirements  for  the  formation  of  persistent  IR  spectral 
holes  (PIRSH's)  that  last  much  longer  than  excited  vibrational  state  lifetimes  at  low 
temperatures?  In  a  general  sense,  the  basic  requirements  for  PIRSH  formation  are  three. 
First,  there  must  be  several  ground  state  configurations  of  the  total  system,  and  the  IR 
energies  from  these  ground  states  must  differ  by  more  than  the  laser  linewidth.  Second, 
there  must  exist  an  IR  pumping  pathway  that  connects  these  ground  state  configurations. 
Finally,  the  relaxation  among  the  ground  states  must  be  slower  than  the  excited  state 
lifetime.  If  all  of  these  conditions  are  met,  persistent  spectral  holes  may  result. 

To  date  persistent  IR  spectral  holes  have  been  found  in  solids  with  three  different 
kinds  of  bonding:  van  der  Waals[40],  ionic[41-43]  and  covalent[44,45]  These  new  high 
resolution  measurements  are  attracting  widespread  interest  not  only  because  they 
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complement  earlier  results  obtained  with  other  spectroscopic  techniques  on  defects  in  solids 
but  also  because  the  source  of  the  persistence  is  often  not  understood.  A  brief  description 
of  persistent  spectral  hole  burning  results  in  the  vibrational  spectrum  of  CN'  and  N02'  in 
alkali  halide  crystals  and  HS  in  As2S3  is  given  below. 

3.2.2  CN-  in  alkali  halide  crystals 

Figure  10  shows  the  low  temperature  absorption  coefficient  in  the  stretch  mode 
frequency  region  for  a  KBr  crystal  doped  with  0.05  mole  %  KCN.  The  vertical  dashed 
line  in  the  figure  divides  the  spectrum  up  into  two  fiequency  intervals.  The  absorption 
lines  which  occur  at  frequencies  larger  than  2100  cm'1  are  associated  with  different 
isotopes  of  the  linear  NCO\  a  common  impurity  in  the  KCN  dopant.  The  linewidths  of 
these  transitions  are  all  about  0.05  cm'1  (FWHM).  None  of  these  features  show  spectral 
hole  burning. 

The  absoiption  lines  at  frequencies  below  2100  cm-i  in  Fig.  10a  are  all  associated 
with  the  CN'  stretch  vibration  in  some  way  or  other  and  the  complexity  of  this  spectrum  is 
quite  surprising.  The  fine  structure  near  the  main  02C14N')  absorption  line  can  be  seen 
more  clearly  in  Fig.  10b  where  the  spectrum  of  a  lower  concentration  sample  (0.01  mole 
%  KCN  )  is  displayed.  This  spectrum  is  characterized  by  one  strong  central  line  with  a 
large  linewidth  (~  0.2  cm'1  FWHM)  surrounded  by  lines  which  have  much  narrower 
linewidths  (~  0.04  cm'1  FWHM).  In  the  figure  these  satellite  lines  are  labelled  ctj,  a^,  ylt 
y2,  Yj  and  Y4-  Our  hole  burning  measurements[42]  indicate  that  both  the  central  line  and 
related  broad  structure  should  be  attributed  to  tunneling  CN'  and  the  ultrasharp  lines,  to 
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CN'  molecules  frozen  in  direction  by  other  nearby  impurities  in  the  crystal. 

No  hole  burning  has  been  detected  in  the  main  line  shown  in  Fig.  10b  but  the  sharp 
lines  labelled  ot]  ,  02  ,  Yi ,  Y2  >  Y3  and  7 4  all  show  persistent  changes  when  probed  with  a 
diode  laser.  Examples  of  hole  burning  in  the  three  absorption  lines  at ,  ,  and  are 

shown  in  Fig.  11.  Figure  11a  shows  a  hole  burned  in  the  dj  line  which  can  made  100  % 
deep  but  which  washes  out  after  about  three  minutes  in  the  dark.  The  lines  y2  ,  Y3  and  y4 
(not  shown)  display  a  similar  erasing  behavior.  Only  a  shallow  hole  can  be  burned  in  the 
(*2  line  [shown  in  Fig.  lib]  which  erases  in  about  30  sec.  Figure  11c  shows  the y{  mode 
before  and  after  two  holes  have  been  burned  in  it.  These  holes  which  can  be  made  100% 
deep  have  the  interesting  property  that  they  are  permanent  as  long  as  the  sample  is 
maintained  at  liquid  helium  temperatures. 
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The  temperature  stability  of  the  holes  burned  in  the  Yi  band  has  been  determined  by 
sweeping  the  spectrum  with  a  low  intensity  probe  beam  while  the  sample  temperature  is 
varied.  Scans  at  four  different  temperatures  are  shown  in  Fig.  12:  trace  (a)  of  a  sample  at 
1.7  K  displays  a  spectrum  which  contains  two  holes,  trace  (b),  the  same  sample  at  15  K, 
trace  (c),  at  25  K  and  trace  (d),  at  35  K.  At  this  highest  temperature  the  spectrum  appears 
nearly  featureless  yet  upon  recooling  to  1.7  K  the  original  spectral  shape  is  recovered. 
Erasing  does  not  become  noticeable  until  the  sample  temperature  is  cycled  up  to  about  48 
K. 

A  systematic  study  of  crystals  double  doped  with  CN'  and  with  different  cation  and 
anion  dopants  demonstrates  that  the  vibrational  center  which  produces  the  Yi  band  is  a  CN' 
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molecule  with  a  near  neighbor  Na+  impurity  replacing  one  of  the  K+  host  ions[42].  Other 


Figure  12. 


complexes  consist  of  the  CN'  molecule  together  with  an  alkali  or  halogen  ion  impurity  at  a 
nearby  lattice  site.  Of  the  fourteen  new  centers  identified,  six  of  them  show  persistent  hole 
burning  when  the  perturbed  CN-  stretch  mode  is  probed  with  a  low  power  diode  laser.  The 
extremely  stable  defect  complex:  CN':  Na+,  which  both  fluoresces  and  hole  bums,  has 
been  investigated  in  some  detail.  By  combining  broadband  Fourier  transform 
interferometry  with  narrow  band  diode  laser  hole  burning  it  has  been  demonstrated  that  the 
particular  persistent  effect  associated  with  the  Yi  band  occurs  when,  during  vibrational  de¬ 
excitation,  the  CN"  molecule  flips  by  180°  between  inequivalent  energy  configurations 
generated  by  the  presence  of  the  nearby  Na+  ion. 

3.2.3  N02* :  the  asymmetric  top  defect 

When  the  nitrite  molecule  N02'  is  doped  into  single  crystal  KI  a  number  of  spectral 
features  appear  in  the  gap  region  of  the  pure  crystal  phonon  spectrum  which  have  remained 
unexplained  for  a  quarter  century.  These  excitations,  which  have  been  studied  by  direct  far 
infrared  (FIR)  absorption[46-48],  infrared  (IR)  sideband  absorption[49],  Raman 
scattering[50-53],  ultraviolet  (UV)  absorption[53,54]  and  emission[55],  have  been 
attributed  to  some  combination  of  translational  modes  and/or  librational  modes  of  the 
molecule  but  no  asignment  is  consistent  with  all  the  available  spectroscopic  data  nor  have 
the  gap  mode  features  ever  been  resolved. 

A  new  solid  state  high  resolution  spectroscopic  technique  has  been  developed  to 
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provide  information  about  the  source  of  the  gap  mode  features  in  the  KI:  N02*  system:  it 
involves  PIRSH  burning  in  the  IR  while  probing  concommitant  changes  in  the  FIR. 
Initially,  tunable  infrared  diode  lasers  are  used  to  burn  persistent  IR  spectral  holes 
(PIRSHs)  in  the  internal  vibrational  spectrum  of  N02"  defects  in  KI  at  low  temperatures 
through  a  reorientational  form  of  hole  burning.  Next,  the  technique  of  PIRSH  burning  in 
the  internal  bending  mode  of  N02‘  is  combined  together  with  high  resolution  FIR  fourier 
transform  interferometry  to  probe  simultaneously  the  KI:N02'  gap  features.  We  find  that 
persistent  polarized  changes  in  the  FIR  spectrum  are  coupled  to  persistent  changes  in  the 
internal  vibrational  mode  spectrum  but  the  production  of  persistent  holes  in  the  IR  does  not 
translate  into  FIR  holes  in  the  external  modes.  The  results  show  that  the  three  FIR 
absorption  lines  in  the  gap  are  consistent  with  localized  translational  motion  of  the  N02* 
ion[43]. 

The  experimental  arrangement  for  probing  the  crystal  simultaneously  with  the  diode 
laser  pump  beam  and  the  FIR  interferometer  probe  beam  is  shown  in  Fig.  (13).  FIR 
transmission  spectra  are  obtained  with  a  Bomem  interferometer  with  a  maximum  apodized 
resolution  of  0.004  cm-1  (120  MHz).  The  interferometer  is  configured  with  a  Hg  arc 
source,  12  mm  mylar  beam  splitter,  and  a  Si  composite  bolometer  at  4.2  K.  The 
interferometer  beam  is  focused  into  a  sample  located  in  a  superfluid  He  immersion  optical 
access  cryostat  (not  shown).  The  laser  polarization  is  rotated  to  any  angle  using  a  quarter 
wave  plate  and  linear  polarizer.  Efficient  illumination  of  the  1.5  x  5.0  mm  sample  face  seen 
by  the  IR  beam  is  accomplished  with  a  ZnSe  cylindrical  lens.  Metal  film  attenuators  are 


Figure  13. 
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used  to  decrease  the  laser  intensity  after  burning  so  that  the  persistent  hole  shape  can  be 
monitored  without  distortion.  Verification  of  proper  IR  sample  illumination  and  monitoring 
of  the  persistent  spectral  hole  is  provided  by  a  HgCdTe  or  a  Ge:Cu  detector. 

Figure  (14)  shows  the  PIRSH  burning  effects  that  can  be  produced  with  a  low 
power  diode  laser  in  the  v2  (A1  symmetry)  infrared  bending  mode  of  N02*  at  804.85  cm*1. 
Curve  (a)  is  the  undistorted  absorption  band  as  measured  in  transmission  by  tuning  the 
attenuated  laser  beam  through  a  few  GHz  about  the  line  center.  To  bum  a  hole,  the  laser 
frequency  is  returned  to  the  center  of  the  absorption  band  and  held  fixed.  At  the  highest 
intensities  used  (25  mW/cm2)  the  absorption  at  the  laser  frequency  decreases  to  zero  in 
about  ten  seconds.  In  the  dark  this  absorption  hole  decays  roughly  exponentially  with  a 
decay  time  of  50  minutes  at  1.5  K.  Holes  burned  with  an  intensity  less  than  0.2  mW/cm2 
cannot  be  burned  100%  deep,  and  no  detectable  hole  is  burned  with  an  intensity  less  than  2 
mW/cm2.  The  hole  shown  in  Fig.  (14),  curve  (b)  is  about  200  MHz  wide  only  because  the 
hole  is  100%  deep.  Note  that  although  antihole  wings  appear  they  are  much  too  small  to 
account  for  the  missing  area  associated  with  the  appearance  of  the  spectral  hole.  Holes  as 
narrow  as  37  MHz  have  been  burned  in  this  band  in  the  low  intensity,  short  bum  time 
limit. 


Figure  14. 


The  PIRSH  burning  process  is  summarized  as  follows:  The  laser  excites  the 
vibrational  bending  mode  of  an  N02*  ion  originally  in  one  of  the  12  allowed  orientations. 
During  de-excitation  an  ion  reorients  with  low  probability  (  ~  10*2)  to  one  of  its  other  1 1 
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positions.  Three  of  these  preserve  the  dipole  axis  (not  necessarily  direction),  and  these 
three  antihole  positions  appear  in  the  same  polarization.  The  other  8  positions  have  their 
permanent  dipole  directions  orthogonal  to  the  starting  direction.  The  long  time  branching 
ratios  for  reorientation  appear  to  be  roughly  equal.  Barriers  to  reorientation  at  1.5  K  allow 
the  molecular  ions  to  remain  in  a  nonthermal  equilibrium  distribution  of  orientations  for  an 
extended  period  of  time  after  the  burning  process.  Since  PIRSH  burning  effects  have 
never  been  investigated  in  the  FIR,  we  now  proceed  to  examine  the  effect  on  the  gap  mode 
spectrum  of  KI:N02_. 

The  low  temperature  absorption  spectrum  of  KLNO2"  in  the  gap  mode  region  at  an 
interferometer  resolution  of  0.015  cm1  is  shown  in  Fig.  (15a).  The  doublet,  which  is 
resolved  here  for  the  first  time,  increases  the  number  of  known  strong  IR  active  N02' 
modes  in  the  gap  from  two  to  three.  In  addition,  each  mode  is  found  to  have  a  weak  low 
frequency  shoulder  consistent  with  an  isotope  shift  due  to  the  7%  abundant  41K  isotope 
in  the  KI  host. 

When  the  internal  bending  mode  is  IR  hole  burned,  the  FIR  spectrum  changes  and 
these  results  are  summarized  in  Figures  (15b)  and  (15c).  In  the  arrangement  shown  in  Fig. 
(13) ,  the  diode  laser  beam  is  unblocked  and  an  IR  hole  is  burned  in  the  v2  mode.  After  30 
minutes  of  burning  at  an  intensity  of  about  0.2  mW/cm2  to  reach  a  steady  state  hole,  the 
FIR  spectrum  is  obtained  with  an  integration  time  of  about  an  hour  at  a  resolution  as  high 
as  0.015  cm1  (450  MHz)  as  determined  from  the  width  of  low  pressure  water  vapor  lines. 
Next  the  diode  is  blocked  and  the  sample  is  warmed  to  10  K  to  erase  the  hole,  cooled  to 
1.5  K  again  and  an  FIR  reference  spectrum  is  obtained.  Ratioing  the  burned  to  the 
unbumed  sample  spectra  reveals  the  absorption  coefficient  changes  for  the  N02‘  molecule 
shown  in  the  figure.  The  changes  produced  in  these  external  modes  do  not  display  sharp 
holes  but  instead  consist  of  line  strength  changes  over  the  entire  FIR  lines.  Erasing  the 
entire  IR  bending  mode  in  each  of  these  experiments  only  increases  the  strength  of  these 
FIR  changes.  In  Figures  (15b)  and  (15c)  the  laser  polarization  is  configured  either  parallel 
or  perpendicular  to  the  interferometer  beam  axis,  hence  a  hole  is  burned  in  v2  either  parallel 
or  perpendicular  to  this  axis.  In  (15b)  the  laser  reorients  some  of  the  permanent  dipoles 
lying  along  the  FIR  axis  into  the  normal  plane.  In  (15c)  the  laser  reorients  some  of  the 
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permanent  dipoles  that  are  normal  to  the  FIR  axis  into  that  direction.  The  observed  changes 
are  a  clear  indication  that  the  transition  dipole  moments  of  these  external  modes  are 
polarized  and  have  a  fixed  relationship  relative  to  the  internal  molecular  axes.  When  the 
permanent  dipoles  are  rotated  from  the  FIR  axis  to  normal  to  the  FIR  axis  as  in  (15b),  the 
71cm*1  mode  is  observed  to  increase.  The  converse  is  also  true,  as  is  shown  in  Fig.  (15c). 
The  sign  and  magnitude  of  the  changes  in  the  71  cm'1  mode  are  consistent  with  a  transition 
dipole  that  is  along  the  permanent  dipole  axis  of  the  molecule,  hence  this  mode  has 
symmetry.  The  higher  frequency  changes  indicate  that  the  polarized  transition  moments  of 
the  79  cm*1  modes  are  perpendicular  to  the  permanent  dipole  axis  (hence  Bj  or  B2 
symmetry). 


In  the  past,  one  of  the  arguments  against  a  local  mode  assignment  has  been  that 
only  two  modes  were  observed  in  the  FIR.  Since  now  three  gap  modes  have  been  resolved 
and  since  the  polarized  PIRSH  burning  results  are  consistent  with  those  expected  for  a 
constrained  Qv  molecule  in  a  cubic  site,  it  is  reasonable  to  assign  these  modes  to  the  three 
expected  nondegenerate  translational  species:  Ax,  Bj,  and  B2  .  All  evidence  on  the 
symmetry  of  the  71  cm-1  gap  mode,  i.e.,  it  alone  appears  in  combination  on  UV  side  band 
spectra  and  has  been  established  to  have  its  transition  dipole  axis  parallel  to  the  permanent 
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dipole  axis  in  this  work,  indicates  that  it  must  transform  with  the  A]  representation. 

One  unanswered  question  has  to  do  with  the  lifetime  of  these  gap  modes.  Our  high 
resolution  FTR-ER  holebuming  measurements  can  be  interpreted  in  either  of  two  ways:  (1) 
The  N02'  gap  modes  are  homogeneously  broadened  and  the  dephasing  or  T2  time  is 
shown  in  Fig.  (15)  or  (2)  These  modes  are  inhomogeneously  broadened  by  local  strains 
with  the  frequency  positions  of  burned  homogeneous  subsets  for  the  internal  bending  mode 
and  external  gap  modes  uncorrelated.  This  interesting  question  could  be  settled  if  a  tunable 
low  intensity  FIR  source  with  a  spectral  width  narrower  than  0. 1  cm1  was  available  with 
which  to  bum  persistent  holes  directly  in  these  gap  modes. 

With  the  combined  PIRSH  burning  plus  FIR  detection  it  has  been  possible  to 
determine  the  isotope  effect  for  the  different  gap  modes.  This  measurement  can  be  used  to 
addressed  a  fundamental  question  in  lattice  dynamics  theory  which  is  now  briefly 
described.  The  vibrational  properties  of  N02‘  in  KI  have  long  attracted  interest  because  of 
the  apparent  conflicting  signature  of  this  substitutional  defect  when  measured  with  different 
spectroscopic  probes.  One  striking  property  is  the  orthogonality  between  direct  far  infrared 
absorption  (FIR)  and  the  Raman  scattering  results  produced  by  this  defect  for  impurity 
modes  in  the  gap  between  the  acoustic  and  optic  phonon  branches  of  the  host 
crystal.  [5 1,53]  To  account  for  this  puzzling  result  Evans  and  Fitchen[53]  noted  that  since 
parity  appears  to  be  a  good  quantum  number  the  C2v  defect  must  not  distort  the  cubic 
symmetry  of  the  I"  cavity.  On  the  other  hand,  parity  is  not  conserved  in  the  model  of 
Rebane  et  al.[51]  where  the  detailed  lattice  dynamical  calculations  of  Zavt[56]  are  used  to 
characterize  the  translational  and  librational  motion  of  the  molecule  in  the  lattice.  Their  fit 
has  been  made  with  anisotropic  force  constants  at  the  defect  site  so  that  the  three  fold 
degenerate  translational  mode  associated  with  an  isotropic  defect  is  split.  Two  translational 
gap  modes  bracket  a  third  mode  which  has  a  mixed  translational  and  librational  character. 
In  their  model  it  is  this  librational  character  which  causes  the  different  spectroscopic  results 
since  it  makes  the  intermediate  mode  strongly  Raman  active  and  only  weakly  FIR 
active[51].  The  possibility  that  this  spectrum  is  mainly  translational[46,47]  or 
librational[49]  in  character  or  even  that  the  Raman  and  FIR  results  obey  the  parity  selection 
rule  just  as  for  a  point  defect[53]  in  a  cubic  lattice  is  still  unresolved. 
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We  have  made  the  first  identification  of  the  isotope  effect  on  the  FIR  gap  mode 
spectrum  by  PIRSH  burning  and  establish  it  as  an  important  technique  for  assigning 
spectral  lines  to  a  particular  isotopic  species  of  N02*.  The  substitution  of  180  for  ,60  and 
15N  for  14N  in  N02*  produces  gap  mode  frequency  shifts  for  two  of  the  three  modes 
consistent  with  translational  mode  behavior.  The  15N  for  14N  substitution  also 
demonstrates  that  the  intermediate  mode  for  the  most  abundant  isotopic  species  does  have 
some  librational  character  even  though  it  does  not  appear  in  the  Raman  spectrum.  This  new 
level  of  FIR  sensitivity  confirms  that  the  two  spectroscopic  techniques  do  probe  orthogonal 
aspects  of  the  defect  dynamics. 

Figure  (16)  shows  the  impurity  induced  spectrum  obtained  for  a  KI  crystal  doped 
with  0.1  mole  %  KNO3  with  40%  l80  enrichment.  The  lines  near  73,  78  and  87  cm1 
appear  to  be  isotopically  shifted  NOf  modes.  The  broad  features  at  92  and  94  cm1  do  not 
show  an  oxygen  isotope  shift.  PIRSH  burning  on  the  v2  bending  mode  of  the  (14,16,18) 
species  with  laser  polarization  parallel  and  perpendicular  to  the  interferometer  beam  axis 
produces  the  FIR  spectal  changes  shown  in  Fig.  (16b)  and  (16c).  Traces  (16d)  and  (16e) 
show  the  gap  mode  changes  produced  by  burning  the  (14,18,18)  IR  bending  mode  with 
the  same  two  laser  polarizations.  Our  investigation  indicates  that  only  spectral  lines  within 
those  grouped  near  71  and  79  cm1  in  Fig.  (16)  are  induced  by  the  different  isotopes  of  the 
NC^'  molecule.  Without  the  PIRSH  labelling  technique,  unambiguous  mode  identification 
would  be  extremly  difficult  and  tedious  in  this  mixed  mode  spectrum. 
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3.2.4  HS  in  As2S3  glass 

All  persistent  spectroscopic  effects  observed  to  date  in  glassy  hosts  have  involved 
electronic  excitation  of  the  defect  center,  with  the  hole  burning  attributed  either  to 
photochemical  modification  of  the  defect  itself,  or  to  rearrangement  of  the  host  atoms 
around  the  defect  Because  of  the  electronic  excitation,  however,  the  observed  behavior  is 
not  necessarily  representative  of  the  glass  near  its  ground  state  configuration.  By  using  a 
vibrational  excitation  in  the  electronic  ground  state  in  the  hole  burning  process,  it  should  be 
possible  to  study  behavior  which  is  more  representative  of  the  pure,  unperturbed  glass. 
Since  persistent  IR  spectral  hole  (PIRSH)  burning  through  defect  reorientation  has  been 
observed  for  vibrational  modes  of  some  defects  in  ionic  crystals  and  van  der  Waals 
matrices[57]  a  search  was  initiated  for  a  similar  effect  in  glasses. 

Here  the  first  observation  of  PIRSH  burning  of  an  impurity  vibrational  mode  in  a 
glassy  host  is  outlined.  The  defect-host  system  is  the  SH  stretching  mode  in  hydrogen- 
doped  As2S3  glass  and  the  results  demonstrate  that  spectral  persistence  can  occur  in 
covalently  bonded  glasses  even  when  the  system  remains  in  the  electronic  ground  state.  In 
addition  to  its  non-electronic  nature,  this  newly  discovered  PIRSH  system  displays  unique 
properties  which  distinguish  it  from  other  systems.  First,  the  large  asymmetry  of  the 
inhomogeneous  broadening  mechanism  for  this  defect  allows  the  connection  between  the 
inhomogeneous  broadening  perturbation  and  the  excited  state  lifetime  to  be  studied. 
Second,  the  spectral  holes  become  narrower  as  they  spontaneously  refill,  which,  to  our 
knowledge,  has  not  previously  been  observed  for  any  system.  Finally,  the  convenient 
timescale  on  which  this  refilling  occurs  permits  it  to  be  monitored  over  a  much  wider 
dynamic  range  than  in  previous  work,  leading  to  a  more  accurate  determination  of  the 
functional  form  of  the  hole  decay,  which  is  found  to  be  well  described  by  a  Gaussian 
distribution  of  tunnelling  parameters. 

Figure  17  shows  the  infrared  absorption  spectra  of  two  samples,  one  doped  with 
hydrogen,  the  other  with  deuterium.  The  hydrogen-doped  sample  displays  an  absorption 
band  peaked  at  2485  cm'1  with  a  full  width  at  half  maximum  (FWHM)  of  74  cm-1  at  1 .5  K. 
Note  that  the  band  width  is  much  larger  than  observed  for  vibrational  defects  in  alkali 
halides.  In  the  deuterated  sample  the  absorption  maximum  is  at  1807  cm1  and  the  width  is 
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51  cm*1.  The  isotope  shift  by  a  factor  of  1.375,  approximately  the  square  root  of  the  ratio 
of  the  masses,  confirms  that  these  absorptions  are  vibrational  in  nature.  The  sharp  line  at 
2323  cm-1  is  the  asymmetric  stretch  mode  of  C02  molecules  disolved  in  the  glass.  PIRSH 
burning  has  recently  been  observed  for  this  defect  as  well  as  for  several  others,  showing 
that  PIRSH  burning  is  a  property  of  the  glass  and  does  not  rely  strongly  on  the  peculiarities 
of  the  specific  defect. 


Figure  17. 


1 600  2000  2400 

Frequency  (cm-1) 


Both  the  SH  and  the  SD  absorption  bands  exhibit  persistent  hole  burning,  and  the 
behavior  of  the  two  isotopes  is  qualitatively  similar.  Holes  have  been  burned  throughout 
the  SH  band  at  temperatures  ranging  from  1.5  K  up  to  47  K.  All  holes  can  be  erased  by 
heating  above  60  K  for  a  few  minutes. 

In  a  typical  amorphous  hole  burning  system,  spectral  diffusion  causes  the  hole 
width  to  increase  with  time  after  burning  has  ceased[58,59].  As  can  be  seen  in  Fig.  18, 
however,  in  the  As2S3:HS  system  the  hole  actually  becomes  narrower  as  it  refills.  The 
simplest  explanation  for  the  hole  narrowing  behavior  is  that  the  excited  state  lifetimes  and 
dephasing  rates,  and  hence  the  homogeneous  widths,  are  not  the  same  for  all  centers 
absorbing  at  a  given  frequency,  and  that  those  centers  with  the  largest  homogeneous 
linewidths  return  most  rapidly  to  the  original  unbumed  configuration.  The  effect  is  seen 
only  for  relatively  long  bum  times  (several  minutes)  or  high  intensities;  the  broad,  quickly 
relaxing  contribution  to  the  hole  takes  longer  to  build  up  than  the  slowly  relaxing  narrow 
contribution.  Since  the  time  evolution  of  the  hole  width  is  the  opposite  of  that  expected  for 
spectral  diffusion,  it  seems  likely  that  spectral  diffusion  does  not  contribute  significantly  to 
the  hole  width.  Hence  the  measured  width  in  the  low  intensity,  short  bum  time  limit 
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reflects  the  true  homogeneous  width,  which  for  Lorentzian  lineshapes  is  half  the  hole 
width. 


Figure  18 


A  striking  feature  of  the  SH  and  SD IR  absorption  bands  is  their  highly  asymmetric 
shape,  with  a  sharp  high  frequency  cutoff  and  long  low  frequency  tail.  This  lineshape  has 
been  observed  to  be  concentration-independent  for  concentrations  ranging  over  almost  three 
orders  of  magnitude,  thus  ruling  out  interactions  between  defects  as  a  factor  determining 
the  lineshape.  At  the  high  frequency  limit  of  these  bands,  the  SH  and  SD  stretch 
frequencies  in  As2S3  are  only  about  2  percent  lower  than  published  gas  phase  values  for 
free  SH  and  SD  radicals[60-62]  (2598.77  cm-1  and  1885.8  cm-1  ,  respectively).  The 
absorption  maxima,  however,  occur  at  a  shift  of  about  4  percent  from  the  gas  phase  values, 
and  the  low  frequency  tails  extend  to  shifts  of  several  times  this  amount.  Since  all 
observed  SH  absorption  in  the  glass  occurs  at  frequencies  lower  than  for  free  molecules, 
and  since  there  exists  a  sharp  high  frequency  cutoff,  it  seems  reasonable  to  interpret  the 
asymmetric  shape  of  the  inhomogeneous  band  as  arising  from  a  perturbation  which  can 
only  act  to  reduce  a  center's  vibrational  frequency.  Absorption  on  the  high  frequency  side 
of  the  band  would  then  be  due  to  relatively  unperturbed  centers,  while  strongly  perturbed 
centers  would  give  rise  to  the  low  frequency  tail  of  the  band. 

If  this  picture  is  correct,  then  an  SH  defect's  frequency  shift  downward  from  the 
high  frequency  cutoff  is  a  measure  of  the  host-defect  coupling  strength  for  that  site.  Since 
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both  the  optical  dephasing  rate  and  the  excited  state  lifetime  depend  on  how  strongly  the 
defect  is  coupled  to  the  host,  one  would  expect  larger  dephasing  rates,  shorter  lifetimes, 
and  hence  greater  hole  widths,  the  lower  the  bum  frequency.  Figure  19  shows  that  this  is 
indeed  the  case,  and  the  effect  is  quite  dramatic.  Plotted  is  the  hole  width  as  a  function  of 
bum  frequency  in  the  short  bum  time  limit  for  holes  burned  at  1.5  K.  For  bum  frequencies 
above  the  absorption  peak,  the  hole  width  is  roughly  constant,  varying  from  0.8  GHz  at 
2525  cm*1  to  0.9  GHz  at  2497  cm-1.  Below  this,  however,  the  hole  width  increases 
sharply  with  decreasing  bum  frequency,  growing  by  an  order  of  magnitude  by  2445  cm1 
and  eventually  exceeding  the  SDL's  0.5  cm-1  continuous  tuning  range. 


In  conclusion,  the  observation  of  PIRSH  burning  for  the  SH  defect  in  As2S3 
shows  that  vibrational  hole  burning  is  not  restricted  to  ionic  and  van  der  Waals  solids,  but 
can  occur  in  covalently  bonded  glasses  as  well.  The  time  dependence  of  the  spontaneous 
hole  relaxation  is  well  described  by  a  Gaussian  distribution  of  tunnelling  parameters,  but  is 
strongly  dependent  upon  the  intensity  used  to  bum  the  hole.  The  temperature  dependence 
of  the  hole  width  suggests  that  the  low  temperature  width  is  dominated  by  the  excited  state 
lifetime.  The  remarkable  increase  in  hole  width  with  decreasing  bum  frequency  is  thus  a 
consequence  of  a  sharp  reduction  of  the  excited  state  lifetime  with  increasing  host-defect 
coupling. 
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3.2.5  Persistent  hole  burning  in  the  far  infrared  (first  attempt) 

Our  successful  demonstration  of  PIRSH  burning  for  N02'  in  alkali  halides  and  the 
observation  of  persistent  changes  in  the  FIR  when  burning  in  the  IR  leads  one  to  ask 
whether  persistent  effects  can  be  produced  by  far  infrared  radiation  itself.  The  fundamental 
experimental  problem  is  the  absence  of  a  tunable  FIR  source.  Although  such  a  source  is 
required  before  such  experiments  can  be  attemped  for  a  variety  of  defect-lattice  systems,  it 
has  been  possible  to  carry  out  such  an  experiment  for  the  gap  modes  produced  by  N02~ 
with  a  fixed  frequency  source. 

In  our  first  attempt  to  try  FIR  persistent  holebuming  experiments  we  collaborated 
with  Professor  V.  Jaccarino  at  the  FIR  Free  Electron  Laser  facility  at  UC  Santa  Barbara  to 
make  use  of  the  tunability  of  their  source.  Although  it  is  true  that  this  laser  can  be  tuned  to 
some  extent  it  has  not  been  possible  to  identify  coincidences  with  the  N02~  gap  modes  in 
KI.  Failing  there  we  attacked  the  problem  at  Cornell  with  a  brute  force  experimental 
method.  The  approach  was  first  to  catalogue  all  the  known  fixed  frequency  gas  laser  lines 
in  the  FIR  and  then  to  measure  the  mode  positions  for  different  isotopic  combinations  of 
N02*  to  identify  any  accidental  coincidences.  Since  the  gap  modes  have  spectral  widths  of 
order  0.1  cm1  a  more  accurate  determination  of  nearby  laser  lines  was  required  than  was 
available  from  the  published  literature. 

The  measurement  of  FIR  laser  frequencies,  good  to  four  significant  figures, 
required  the  construction  of  a  scanning  metal  mesh  Fabiy-Perot  interferometer.  A  side 
view  of  this  instrument  is  shown  in  Fig.  (20).  Metal  meshes  with  spacing  to  wavelength 
such  that  99%  of  the  radiation  is  reflected  are  stretched  over  parallel  metal  rings  to  form  two 
reflectors  in  tandum.  These  free-standing  mesh  reflectors  (A  in  Fig.  20),  made  parallel  by 
adjustment  of  micrometers  (B),  are  translated  relative  to  each  other  in  one  micron  steps  by  a 
computer-controlled  stepper-motor-driven  translation  stage  (C).  Because  of  the  high 
reflectivity  of  the  mesh,  FIR  transmission  through  the  pair  peaks  sharply  whenever  their 
separation  d  is  a  multiple  of  the  half-wavelength  (d  =  n  1/2  where  n  is  an  integer).  This 
instrument  has  been  used  to  measure  wavelengths  to  a  part  in  104  and  FIR  laser  linewidths 
down  to  100  MHz. 
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Figure  20 


These  results  together  with  the  isotope  shifted  gap  modes  with  center  frequencies  identified 
to  five  significant  figures  are  presented  in  Table  2.  Inspection  of  this  Table  shows  that  a 
coincidence  within  a  gap  mode  line  width  still  does  not  occur  accidentally. 

In  order  to  improve  the  coincidence  between  the  fixed  laser  frequency  and  the  gap 
mode  we  have  prepared  alloy  samples.  The  linear  dependence  of  the  mean  lattice  constant 
on  the  alloy  concentration  is  called  Vegard's  law.  Since  the  gap  mode  frequencies  also 
depend  on  the  lattice  constant  of  the  host  crystal,  some  frequency  shift  is  to  be  expected  in 
alloy  crystals.  We  have  found  that  this  "poor  man  tuning"  corresponds  to  about  0.1  cm-1 
for  every  2  mole  %  of  KBr  in  KI.  Similar  effects  have  been  found  for  other  combinations. 
If  the  lattice  constant  of  the  minor  component  is  larger  than  the  host  lattice  of  the  major  one 
then  the  gap  modes  move  to  lower  frequencies  upon  alloying  and  they  go  in  the  opposite 
direction  if  the  converse  is  true.  Since  the  alloying  broadens  the  modes  with  a  width 
comparable  to  the  frequency  shift  there  is  little  value  in  using  alloys  above  6%  so  only  a 
small  shift  can  be  realized  in  practice  using  this  technique. 

No  persistent  FIR  holebuming  has  been  observed  for  the  three  coincidences 
studied,  see  the  three  rows  identified  with  a  (*)  in  Table  2.  The  absence  of  any  effect 
indicates  that  if  FIR  hole  burning  does  occur  the  efficiency  is  so  small  that  less  than  one 
photon  in  104  converts  the  molecule  to  the  new  configuration.  For  comparison,  the  PIRSH 
efficiency  for  NO2'  is  10*2  while  the  weakest  optical  hole  burning  efficiency  reported  for 
any  system  is  10*6. 
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Table  2.  FIR  laser  coincidences  with  gap  mode  frequencies  in  KI:  N02*. 


DEFECT 

HOST 

MODE 

FREQUENCY 

(cm"1) 

FIR  GAS 

CW  PUMTCD 
LINE  (cm’1) 

STATUS 

15n16o16o- 

KI 

70.340 

so2 

70.3 

70.23* 

14N180180' 

KI 

70.386 

14N160160' 

KI 

70.463 

14N160180‘ 

KI 

70.730 

15n160160- 

KI 

70.828 

CH^NH^ 

70.9 

PU 

14N160160‘ 

KI 

70.985 

CH3OD 

71 

U 

i5N160160' 

KI 

78.145 

CH3OD 

78.1 

U 

1 

0 

00 

rH 

0 

00 

z 

rH 

KI 

78.204 

15n16o16o* 

KI 

78.417 

•— * 
■u 

5S 

M 

00 

0 

00 

q 

KI 

78.438 

14N160180‘ 

KI 

78.616 

14n16o16o* 

KI 

79.000 

CH3NH2 

79 

78.86* 

14N160180‘ 

KI 

79.015 

ch3nh2 

79 

78.86* 

14n16o18o* 

KI 

79.112 

14N160160' 

KI 

79.124 

15N160160' 

KI 

79.397 

qH4(OH)2 

79.4 

P,U 

15n16o16o* 

KI 

79.510 

C2H4(OH)2 

79.5 

P,U 

14N160160' 

KI 

79.698 

14N16o160* 

KI 

79.804 

14N160160' 

KBr 

97.85. 

CD30H 

97.500 

G 

Legend  to  status  column: 

G:  gas  bought. 

P:  pulsed  CO2  laser  pumped  line . 

U:  tried,  not  found. 
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3.3  Excitations  in  Anharmonic  Crystals  and  Glasses 
3.3.1  Preliminaries 

In  the  past  the  lattice  dynamics  of  anharmonic  crystals  have  been  described  in  the 
same  manner  as  for  harmonic  crystals.  Since  the  atoms  undergo  small  oscillations  the 
important  term  in  the  potential  is  still  harmonic.  In  the  quasi-harmonic  approximation  the 
essential  difference  for  the  anharmonic  case  is  that  the  effective  harmonic  potential  that  each 
atom  sees  is  temperature  dependent  hence  the  frequencies  of  the  plane  wave  phonons  that 
are  used  to  describe  the  collective  vibrational  modes  of  the  system  are  also  temperature 
dependent  These  excitations  are  often  referred  to  as  renormalized  phonons. 

One  of  the  most  anharmonic  crystals  is  solid  3He.  For  such  quantum  systems 
classical  lattice  dynamics  does  not  apply[63,164];  however,  it  is  still  possible  to  construct 
an  effective  potential  energy[65]  so  that  a  semiquantitative  comparison  between  experiment 
and  theory  is  possible.  A  large  number  of  experiments  have  been  carried  out  on  this 
system  to  test  whether  or  not  renormalized  phonons  provide  a  complete  description  of  the 
vibrational  excitations.  The  current  status  is  that  most  data  can  be  described  with  these 
excitations  but  yet  there  are  some  experimental  results  which  are  not  consistent  with  the 
renormalized  phonon  model.  These  experiments  which  demonstrate  the  inadequacies  of  the 
current  theory  are  mentioned  here. 

The  early  measurements  of  the  low  temperature  specific  heat  of  bcc  3He  showed 
three  temperature  dependent  components:  one  which  varied  as  T3  and  a  second  component 
which  varied  linearly  with  temperature  [66-69]  and  at  high  temperatures  a  third  component 
which  varied  exponentially  with  temperature[66,70].  This  situation  prevailed  for  ten  years. 
When  measurements  were  finally  made  with  large  crystals  which  had  been  cooled  very 
slowly  in  an  open  geometry,  the  linear  term  was  no  longer  observed[71]  but  the  other  two 
remained.  These  last  results  appeared  to  indicate  that  the  source  of  the  linear  term  observed 
by  previous  workers  was  extrinsic  to  solid  helium.  This  point  will  be  considered  again 
later. 

The  extra  high  temperature  contribution  to  the  crystal  specific  heat  of  b.c.c.  3He 
was  determined  to  be  of  the  activated  form,  namely, 

Q/R  =  A(ftT)2exp(- ftT) ,  (13) 
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where  f  is  the  activation  energy  for  a  particular  molar  volume.  This  form  for  the  excess 
specific  heat  contribution  does  not  uniquely  identify  the  mechanism  since  it  is  consistent 
with  an  Einstein  oscillator  term  (A  =  3),  a  two  level  Schottky  anomaly  (A  =  gi/go)  where 
gt  and  go  are  the  degeneracy  of  the  excited  state  and  the  ground  state,  respectively,  vacancy 
formation  (A  =  1),  vacancies  and  interstitials  individually  (A  =  4),  or  vacancy-interstitial 
pairs(A  =  3.5)[67,71].  The  original  data  was  fit  with  A  =2  but  the  most  recent  and  most 
accurate  heat  capacity  measurements  show  that  the  contribution  is  best  fit  with  A  ®  1,  that  A 
is  temperature  dependent,  and  that  f  ~  0D2,  where  0D  is  the  Debye  temperature  for  a 
particular  molar  volume  at  0  K  [71].  The  idea  that  this  effect  is  associated  with  vacancy 
formation  (A  =  1)  has  been  pursued  by  a  number  of  workers  [71-73]. 

A  direct  determination  of  an  activated  volume  change  in  3He  has  been  obtained  from 
measurements  of  the  temperature  dependence  of  the  x-ray  lattice  parameter  in  a  constant 
volume  cell[74,  75].  The  thermal  vacancy  concentration  is  equal  to  the  change  in  the 
macroscopic  volume  of  a  cell  minus  the  microscopic  x-ray  volume.  With  this  number 
determined,  it  has  been  shown  that  the  localized  vacancy  assumption  gives  too  large  a 
specific  heat  contribution  to  be  consistent  with  the  measured  data.  Together  the  specific 
heat  and  x-ray  data  demonstrate  that  nonlocalized  vacancies  may  occur. 

Glasses  also  show  anomalous  temperature  dependencies  in  the  thermodynamic 
properties.  In  addition,  as  a  function  of  temperature  the  low  temperature  specific  heat  of 
glasses  shows  a  universal  form[76,77].  At  low  temperatures  the  specific  heat  varies 
linearly  with  temperature  while  for  higher  temperatures  the  T3  dependence  is  observed  but 
the  coefficient  of  this  Debye  specific  heat  term  is  10%  to  100%  larger  than  the  value 
calculated  from  the  measured  sound  velocities  in  the  same  material[77,78]. 

An  early  explanation  of  the  linear  term  made  use  of  tunneling  of  atoms  or  groups  of 
atoms  in  the  disordered  solid.  An  ad  hoc  assumption  that  the  density  of  tunneling  states 
was  a  constant  independent  of  energy  was  required  to  explain  the  observed  temperature 
dependence[79,80].  Although  the  tunneling  model  is  still  used  to  describe  the  low 
temperature  experimental  data,  during  the  past  18  years  a  convincing  argument  which 
produces  both  of  the  characteristic  thermal  properties,  namely  the  linear  term  and  the 
enhanced  coefficient  of  the  Debye  term,  has  not  been  forthcoming. 
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3.3.2  Anharmonic  localized  modes  in  crystals  and  glasses 

During  this  grant  period  the  PI  in  collaboration  with  Professor  S.  Takeno  from 
Japan  have  come  up  with  a  new  model  to  describe  the  dynamics  of  anharmonic  crystals  and 
glasses.  The  proposal  was  described  in  a  series  of  four  concise  and  fairly  mathematical 
papers[81-84].  Here  the  model  is  described  in  a  qualitative  manner  to  bring  out  the  key 
features  which  can  account  for  some  of  the  unexplained  experimental  data  described  in  the 
previous  section. 

To  visualize  how  localized  modes  can  appear  in  a  perfect  but  anharmonic  lattice  the 
one  dimensional  harmonic  lattice  is  considered  first.  The  normal  modes  are  described  by 
an  orthogonal  set  of  plane  homogeneous  waves  which  are  spread  through  a  band  of 
frequencies.  An  excitation  localized  at  a  particular  site  can  be  constructed  from  a  linear 
superposition  of  these  normal  modes  and  this  localized  wave  packet  centered  around  a 
particular  frequency  in  the  band  moves  through  the  crystal  in  the  usual  way. 

For  a  one  dimensional  anharmonic  lattice  in  the  limit  of  small  oscillations  the  plane 
wave  spectrum  looks  the  same  as  for  the  harmonic  case.  What  happens  if  one  constructs  a 
wave  packet  to  represent  a  localized  vibration  in  the  chain?  To  be  specific  let  us  assume 
that  the  anharmonicity  is  quartic  with  a  positive  coefficient.  As  the  wave  packet  is 
constructed  the  amplitude  of  the  atom  increases  and  since  the  potential  is  not  harmonic  so 
does  its  frequency  of  vibration.  Now  as  one  localizes  an  excitation  more  and  more  around 
a  particular  atom  the  amplitude  of  vibration  gets  larger  and  the  center  frequency  moves  up 
through  the  band  of  allowed  frequencies  until  it  pops  out  of  the  top  of  the  spectrum.  Given 
sufficient  anharmonicity  the  end  result  is  a  true  localized  mode  centered  at  the  atom  in 
question.  The  anharmonic  system  has  evolved  in  such  a  way  that  both  homogeneous  and 
inhomogeneous  waves  exist 

The  eigenvector  of  the  highest  frequency  mode  for  three  specific  cases  arc  shown  in 
Fig.  (21).  Case  (c)  depicts  the  highest  frequency  (tom)  plane  wave  mode  of  the  harmonic 
linear  chain.  With  increasing  anharmonicity  a  localized  solution  exists,  (b),  which  is 
separated  off  from  the  plane  wave  spectrum  so  that  there  are  now  only  N-l  plane  wave 
modes  and  one  localized  mode  with  a)„.>ci)m.  For  this  case  there  are  roughly  Nd  atoms 
involved  in  the  localized  vibration  and  since  Nd  «  N  the  amplitude  of  each  atom  in  this 
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mode  is  much  larger  than  for  the  plane  wave  case.  In  the  high  frequency  limit,  £0r  »0)m, 
the  local  mode  eigenvector  takes  on  the  particularly  simple  form  shown  in  case  (a)  with  the 
excitation  completely  localized  on  the  central  atom  and  its  two  neighbors. 
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A  schematic  representation  of  the  frequencies  associated  with  the  plane  wave 
spectrum  located  between  0  and  (0m  and  the  anharmonic  localized  mode  at  co:.  is  shown  in 
Fig.  (22b).  The  relative  frequency  of  the  localized  mode  to  com  depends  on  the  relative 
magnitudes  of  the  anharmonic  and  harmonic  potential  terms.  Also  shown  in  the  picture  is  a 
dashed  line  which  represents  the  position  of  the  local  mode  third  harmonic.  In  principle, 
the  anharmonic  potential  produces  a  response  at  this  frequency,  but  because  it  is  so  much 
larger  than  the  plane  wave  and  the  local  mode  frequencie,  it  is  assumed  that  the  lattice 
cannot  respond  and  so  it  is  ignored  in  the  dynamical  analysis.  This  assumption  is  usually 
referred  to  as  the  rotating  wave  approximation. 

Somewhat  harder  to  visualize  is  the  resonant  mode  which  occurs  in  the  case  of  the 
soft  anharmonic  lattice.  Once  again  a  localized  excitation  results,  but  now  at  low 
frequencies  the  central  atom  vibrates  against  the  rest  of  the  crystal  lattice.  The  characteristic 
frequencies  in  this  problem  are  shown  in  Fig.  (22a).  The  resonant  mode  is  far  down  in  the 
acoustic  spectrum  so  0)r « com.  The  corresponding  eigenvector  is  shown  on  the  right. 
Note  that  for  this  soft  potential  case  a  response  at  the  third  harmonic  is  possible  since  it 
occurs  at  an  energy  in  the  plane  wave  spectrum.  This  additional  channel  shortens  the 
lifetime  of  the  resonant  mode  and  makes  it  less  well  defined  than  for  the  local  mode. 
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Figure  22. 


With  the  fonn  of  the  solutions  identified  for  the  two  kinds  of  potentials  the  next  step 
is  to  introduce  temperature  into  the  problem.  Consider  the  dynamics  of  a  quantum 
mechanical  one  dimensional  harmonic  lattice  of  N  atoms  in  thermal  equilibrium  with  a  heat 
reservoir.  At  absolute  zero  temperature  each  mode  in  the  bounded  plane  wave  spectrum 
has  a  zero  point  energy  of  fico/2.  Because  all  atoms  are  involved  in  each  mode  the  root 
mean  square  amplitude  of  a  given  atom  is  N‘1/2  times  the  value  it  would  have  if  only  this 
one  atom  was  involved  in  this  particular  mode.  At  low  temperature  the  lowest  energy 
modes  of  the  system  within  kT  of  the  zero  point  states  become  thermally  excited  giving  rise 
to  phonons.  Most  of  the  higher  frequency  modes  remain  in  their  zero  point  ground  states 
but  in  any  event  all  excitations  are  still  described  by  plane  waves.. 

Now  consider  what  happens  in  the  case  of  the  anharmonic  crystal.  At  absolute  zero 
temperature  the  plane  wave  description  applies  just  as  for  the  harmonic  case  above. 
However,  the  situation  is  quite  different  at  finite  temperature  since  as  shown  above  both 
homogeneous  and  inhomogeneous  wave  solutions  can  occur.  Given  kT  of  energy  per 
atom  the  system  can  either  distribute  it  over  the  plane  wave  spectrum  producing  phonons  in 
the  usual  way,  or  use  some  of  the  energy  to  create  high  frequency  localized  modes  in  their 
zero  point  motion  ground  states  (for  the  case  of  positive  quartic  anharmonicity)  or  low 
frequency  resonant  modes  (for  the  case  of  negative  quartic  anharmonicity).  These  localized 
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excitations  will  be  referred  to  as  "anons"  in  the  rest  of  the  discussion.  Since  it  takes  energy 
to  create  an  anon  in  its  zero  point  state  this  event  would  be  extremely  unlikely  if  only 
energy  was  involved.  But  the  creation  of  a  localized  defect  in  a  perfect  lattice  also  produces 
a  large  configurational  or  entropy  of  mixing  term  since  the  defect  can  appear  at  any  one  of 
the  N  lattice  sites,  hence  to  determine  their  existence  it  is  important  to  keep  track  of  the  free 
energy  of  the  system.  For  the  one  dimensional  lattice  considered  here  the  entropy  term  is 
still  not  large  enough  to  make  localized  defects  occur. 

Only  in  three  dimensions  does  the  entropy  become  large  enough  to  change  the  sign 
of  the  free  energy  of  the  lattice  so  that  localized  states  automatically  appear  at  finite 
temperature.  The  argument  is  similar  to  that  used  to  account  for  vacancies  in  crystals 
however  the  energy  scale  is  much  smaller  now  since  the  atom  is  not  required  to  move  out 
of  a  given  lattice  site,  instead  the  atom  remains  at  its  equilibrium  position  and  a  localized 
strain  field  accompanies  the  creation  of  each  anon.  The  number  of  anons  in  the  crystal  at 
any  temperature  T  obeys  an  activation  law.  We  have  proposed  that  the  activation  term 
found  in  the  experimental  specific  heat  of  solid  3He  may  stem  from  the  production  of  these 
anon  modes[82]. 

So  far  the  discussion  has  been  limited  to  stationary  anons  in  crystalline  lattices. 
Because  all  lattice  sites  are  equivalent  it  is  expected  that  these  modes  could  translate.  How 
easily  they  move  is  still  an  unanswered  question.  The  answer  must  depend  to  some  extent 
on  the  incommensurability  of  the  length  scale  of  the  shape  function  of  the  anon  with  respect 
to  the  lattice  spacing. 

If  anons  do  not  translate  but  in  fact  diffuse  through  the  lattice,  then  this  motion 
produces  a  contribution  to  the  specific  heat  which  is  linear  in  temperature[84,85]. 
However,  because  of  the  activation  energy  required  to  produce  these  modes,  the  number 
density  is  too  small  to  detect  at  low  temperatures  in  ordinary  dielectric  crystals. 

The  situation  is  very  different  in  glasses  since  disorder  is  frozen  in  upon  cooling 
through  the  glass  transition  temperature.  If  it  is  assumed  that  the  structural  disorder  makes 
local  volume  fluctuations  in  the  solid  which  are  compatible  with  the  requirements  of  the 
anon  configurations,  then  the  maximum  number  of  these  modes  become  available  at  all 
temperatures  below  the  glass  transition  temperature.  In  contrast  with  the  anharmonic  single 
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crystal  where,  because  of  the  activation  energy  associated  with  the  localized  configurations, 
the  low  temperature  modes  are  essentially  all  phonons,  the  disordered  solid  contains  both 
phonons  and  anons.  The  phonons  give  a  T3  contribution  to  the  specific  heat  of  glass  and 
the  diffusing  anons,  a  linear  term. 

The  anons  also  explain  the  fact  that  the  T3  contribution  is  enhanced  over  what  is 
calculated  from  the  velocity  of  sound  measurements.  Assume  that  the  anons  are  associated 
with  resonant  modes.  In  the  limit  that  the  density  of  resonant  modes  is  small,  the  total  low 
temperature  specific  heat  of  the  disordered  solid  is  assumed  to  consist  of  a  sum  of  linearly 
independent  contributions  in  analogy  with  the  defect  resonant  mode  case  in  crystals  so 
CV(T)  =  (n  -  na)cD  +  n^  +  nacw  (14) 

where  the  first  and  last  terms  describe  the  contribution  from  Debye,  and  stationary  resonant 
mode  vibrational  degrees  of  freedom  which  so  far  have  been  ignored.  Note  that  since  the 
the  resonant  mode  is  stationary  it  does  not  contribute  in  a  propagating  pulse  experiment 
such  as  is  carried  out  in  a  sound  velocity  measurement  For  n  »  na  Eq.  (14)  simplifies  to 
CV(T)  *  neD  +  n.Cjd  +  nacav  (15) 

To  estimate  nacav  in  Eq.  (15)  requires  knowledge  of  the  density  of  vibrational  states 
below  the  resonant  mode  center  frequency.  Although  this  quantity  is  not  yet  available  for 
anharmonic  resonant  modes,  the  eigenvalue  equation  for  the  anharmonic  mode  is  similar  to 
that  obtained  for  a  harmonic  resonant  mode  produced  by  a  force  constant  defect.  The  mean 
square  amplitude  of  the  atom  located  at  site  0  versus  frequency  is  shown  for  the  force 
constant  defect  in  Fig.  (23).  If  this  analogy  carries  over  to  the  calculation  of  the  low 
frequency  density  of  states,  which  we  now  assume,  then  for  the  limit  where  to<  o\,  the 
extra  low  frequency  contribution  to  the  vibrational  density  of  states  varies  as  co2  just  as  for 
the  plane  wave  spectrum.  Since  the  form  is  quadratic  in  frequency,  this  vibrational 
component  will  provide  a  T3  contribution  to  the  specific  heat.  The  trade  off  between  the 
increased  density  of  states  at  frequencies  below  the  resonance  and  the  relatively  small 
number  density  of  resonant  mode  centers  is  such  that  a  significant  contribution  to  the  T3 
specific  heat  results.  We  find  that  C^/Cq  =  0.25,  a  constant  in  the  low  temperature  limit. 

The  anon  model  predicts  two  interrelated  signatures  in  the  low  temperature  specific 
heat  of  glasses:  a  T3  term  which  stems  from  an  enhanced  density  of  vibrational  modes  in 


46 


<b) 


—  K, 


-2 


-I  0  I 


2 


Figure  23. 


the  neighborhood  of  <or  and  a  linear  term  which  appears  because  these  resonant  mode 
centers  can  diffuse.  Both  of  these  features  have  been  found  to  be  universal  properties  of 
glasses. 

Because  of  diffusion  the  anon  center  of  mass  motion  cannot  contribute  to  the 
propagation  of  an  energy  pulse  the  way  a  plane  wave  mode  can.  In  addition,  its  localized 
vibrational  degree  of  freedom  corresponds  to  a  standing  wave  so  it  cannot  contribute  either. 
These  characteristics  are  consistent  with  the  idea  that  the  anon  contribution  will  remain 
partially  hidden  in  any  experiment  which  probes  the  propagation  of  waves  through  the  solid 
such  as  a  velocity  of  sound  measurement 
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